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摘要 


[现代 分 形 理论 中 存在 有 许多 问题 尚 待 解决 ， 尤 其 关于 分 维 数 的 计算 、 分 形 反 
问题 及 奇异 性 分 析 等 方面 的 研究 为 重点 。 小 波 在 处 理 非 平稳 性 和 局 域 分 析 方面 有 
着 独特 的 优越 性 ， 可 被 用 来 提取 有 关 分 形 尺度 性 质 的 细微 信息 本文 以 小 波 理论 
为 工具 围绕 分 形 中 若干 问题 进行 了 探讨 ， 主 要 提出 并 解决 了 如 下 问题 ; 

1， 分 析 了 分 形 曲线 维 数 的 传统 计算 方法 及 其 缺陷 ， 提 出 了 多 分 维 分 形 曲 线 
的 模式。 基于 小 波 阔 值 去 噪 方法 ， 给 出 曲线 上 点 的 局 部 分 形 维 数 的 计算 公式 及 算 
法 。 仿 真 计算 结果 及 图 示 均 表明 该 法 是 有 效 的 。 

2， 基 于 小 波 变换 在 零点 具有 的 的 保 相似 形 特 征 ， 在 已 知 分 形 插值 函数 图 象 
情况 下 ， 从 有 限 数目 个 零点 出 发 ， 给 出 寻找 一 组 仿 射 压缩 变换 使 该 函数 图 象 成 为 
未 知 迭 代 函 数 系统 吸引 子 的 一 种 构造 方法 。 

3， 提 出 了 分 形 时 变 维 数 的 概 含 ， 将 其 引入 两 类 具有 局 部 自 相似 性 的 有 偏 随 
机 过 程 中 。 借 助 小 波 基 固有 的 尺度 特性 很 适合 分 析 非 平稳 性 和 局 部 自 相似 性 这 一 
特点 ， 给 出 灰 斯 特 指数 的 小 波 估算 公式 及 算法 ， 证 明了 小 波 估算 信和 真实 值 是 相 
容 的 ， 并 用 仿真 结果 加 以 验证 ， 同 时 对 上 海 股票 市 场 周 指数 变化 进行 了 实证 分 析 。 

4. 借助 小 波 变 换 对 自 相似 过 程 所 起 的 白化 作用 和 最 小 二 乘法 ， 建 立 了 无 须 
知道 V, 过 程 中 参数 分 布 的 谱 参 数 估计 方法 ， 

5， 提 出 在 最 优 尺度 上 跟踪 小 波 模 极 大 值 线 的 参数 方法 ， 为 研究 多 分 形 谱 奇 
异性 提供 了 一 种 途径 。 探 讨 了 分 形 测度 小 波 变 换 极 大 值 线 的 拓扑 分 偏 情况 ， 并 提 
出 了 一 种 基于 小 波 变换 模 极 大 值 的 多 重 分 形 谱 估算 方法 。 
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Abstract 


There are many problems to be solved in modern fractal theory, especially focuses on 
studies of computation on fractal dimension,inverse problem and sigularities. The scale 
property of wavelet is well suited for analyzing non-stationarity and localization and 
used to extract microscopic information about the scaling properties of fractals. In this 
dissertation, based on wavelet theory,the following probloms are discussed: 

1. Traditional computional methods and limitation of dimension of fractal curve 
being analyzed ,A model of fractal curve having multidimension is proposed. Based on 
noise reduction for wavelet softhresholding, computional formula and algorithm about 
point-wise local dimension on the curve are given. Results from simulations show the 
validity of he algorithm. 

2. Given a graph of a fractal interpolation function which is the attractor of an 
unknown IFS with affine constration maps, the maps are found based on the self- 
similarity of the zero-crossing points of wavelet transform. The effectiveness of method 
is shown in an example. 

3. The concept of time-varying dimension are proposed and introduced in two kinds 
of biased stochastic process with locally self-similarity. The estimation formula and 
algorithm of Hurst Index are based on Daubechies wavelet analysis of samples dada . 
Simulation result indicates effectiveness of estimation method and proves that 
estimation value is a consistent result of true value. Week-index analysis of Shanghai 


stock market is taken as a real example by time-varying dimension. 
4. Based on correlation-decaying charecteristics of discrete wavelet coefficients of 


し processes and the least-square algorithm „estimation for spectral parameter of / 


processes is presented without knowing the distribution of the parameter. 

5. A parameterization method for tracing wavelet maxima lines at the finest scale 
continuously is presented, by means of which the sigularity of multifractal spectrum can 
be analyzed well. Topologic bifurcation situation of maxima lines is also discussed. In 
the end ,an estimation for multifractal spectrum is proposed in terms of wavelet module 


maxima lines. 
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第 一 章 绪论 


第 一 章 绪论 
1.1 分 形 理论 的 发 展 及 应 用 


分 形 理论 是 非 线性 科学 研究 中 一 个 十 分 活跃 的 分 支 ， 它 的 研究 对 象 是 目 然 界 非 
线性 科学 中 出 现 的 不 光滑 和 不 规则 的 几何 体 。 昌 然 分 形 理论 在 本 世纪 70 年 代 才 提 出 
来 ， 但 经 过 十 几 年 的 发 展 ， 它 已 广泛 应 用 到 自然 科学 和 社会 科学 的 几乎 所 有 领域 ， 
成 为 当今 国际 上 许多 学 科 的 前 沿 研究 课题 之 一 。 


1. 1. 1 分 形 理论 发 展 的 历史 进程 


第 一 阶段 (1875 年 ---1925 F) ， 人 们 已 经 认识 到 几 种 典型 的 分 形 集 ， 并 且 力 负 
对 这 些 集合 与 经 典 几 何 的 差别 进行 描述 分 类 和 刻画 。19 HARA BA ERE DS 
连续 与 可 微 的 曲线 ， 但 是 普遍 认为 连续 而 不 可 微 的 情形 是 极其 少 的 ， 并 在 理论 研究 
中 应 该 排除 这 类 “怪物 ”， 而 且 更 进一步 认为 一 条 连续 曲线 上 不 可 微 点 应 当 是 很 少 
的 。 维尔 斯 特 拉 斯 (Weierstrass) 在 1872 年 证 明了 一 种 连续 函数 在 任意 一 点 均 不 具有 
有 限 或 无 限 导数 。 这 一 结果 在 当时 曾 引 起 极 大 震动 ， 但 人 们 还 是 认为 维尔 斯 特 拉 斯 
型 的 函数 是 极端 “病态 ”的 例子 。 即使 如 此 ， 人 们 依然 从 不 同 的 方面 推广 了 上 述 馈 
数 ， 并 对 其 奇异 性 作 了 深入 的 研究 ， 获 得 了 极为 丰富 的 结果 。 

康 托 尔 (Cantor) 在 1872 年 引入 了 一 类 全 不 连通 的 紧 集 ， 后 被 称 为 康 托 尔 三 分 集 
趾 。 其 构造 过 程 是 这 样 的 对 于 单位 直线 段 /, ， 去 掉 中 间 三 分 之 一 的 线段 后 所 得 图 
形 为 1 。 然后 对 7 的 每 个 线段 都 进行 同样 的 操作 过 程 。 依次 类 推 ， 便 可 以 得 到 一 个 
曲线 序列 人 }。 当 趋 于 无 穷 大 时 ， 曲线 序列 世上} 趋向 于 一 个 极限 曲线 ， 则 就 是 
康 托 三 分 集 。 这 类 集合 在 三 角 级 数 的 唯一 性 这 样 重要 问题 的 证 明 中 发 挥 了 非常 的 作 
用 。 

1890 年 皮 亚 诺 (Peano) 构 造 出 了 另 一 类 能 够 填充 平面 的 曲线 ， 从 而 引起 人 们 开始 
对 以 入 长度 和 面积 概念 的 反思 凶 。 皮 亚 诺 曲 线 及 其 它 的 例子 导致 了 后 来 拓扑 维 数 概 
念 的 引入 。 

冯 。 科 赫 (Von Koch) 于 1904 年 采用 初等 方法 构造 出 了 现今 被 称 为 汉 。 科 鱼 曲 
线 的 处 处 不 可 微 的 连续 曲线 ， 并 且 讨论 了 曲线 的 性 质 。 这 种 曲线 构造 方面 的 简单 
性 ， 一 改 人 们 一 直 认 为 达 续 不 可 微 曲线 构造 一 定 非常 复杂 的 观念 。 
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由 于 这 一 时 期 各 类 复杂 集合 的 出 现 ， 人 们 重新 认识 长 度 、 面 积 等 这 些 数学 基 概 
念 。 为 了 测量 这 些 集合 ， 在 1901 年 闵可夫 斯 基 (Minkowski) 引 入 了 闵可夫 斯 基 容 上 度 
概念 ， 进 一 步 的 结果 是 ， 豪 斯 道夫 (Hausdorfp 于 1919 年 引入 了 豪 斯 道夫 测度 和 豪 斯 
道夫 维 数 叫 。 这 些 概念 的 产生 表明 ， 要 测量 一 个 几何 对 象 ， 人 们 必须 依赖 于 测量 方 
式 和 測量 所 米 用 的 尺度 。 

在 二 十 世纪 的 初期 ， 一 类 典型 的 随机 分 形 集 一 Brown 运动 各 5 的 已 受到 物理 学 家 
们 的 重视 。1913 年 柏 瑞 (Perrin) 对 粒子 Brown 运动 的 轨迹 图 进行 了 深入 细致 的 研究 ， 
得 出 了 Brown 运动 作为 运动 曲线 不 具有 导数 的 结果 。 他 的 这 一 结论 引起 了 维 纳 
(Wiener) 的 注意 ， 也 直接 促使 了 他 在 1920 年 左右 建立 了 Brown 运动 的 概率 模型 ， 
维 纳 此 时 采用 “混沌 ”一 词 来 表达 这 种 极端 的 混乱 形式 。 

曼 德尔 布 罗 特 (Mandelbrot) 中 在 回顾 这 一 时 期 的 分 形 几 何 历史 时 指出 ， 分 形 几 何 
以 下 面 两 种 选择 方式 作为 其 特征 ， 一 是 在 自然 界 的 混沌 中 选择 问题 ， 二 是 在 数学 中 
选择 问题 。 这 两 种 选择 逐渐 融合 ， 从 而 在 无 序 混沌 和 过 分 有 序 的 欧 氏 几何 之 闻 产 生 
出 了 一 个 具有 分 形 序 的 新 领域 。 

第 二 阶段 (1926 年 一 1975 年 ) ， 在 这 半 个 世纪 里 ， 人 们 对 分 形 集 的 维 数 的 研究 
获得 了 丰硕 的 成 果 b5t7zinslte]， 不 仅 逐 渐 形 成 了 理论 ， 而 且 研 究 的 范围 扩大 到 数学 许 
多 分 支 中 去 。 以 贝 西 可 维 奇 (Besicovitch) 为 代表 的 一 批 学 者 研究 了 曲线 的 维 数 、 分 形 
集 的 结构 及 局 部 性 质 以 及 在 数论 、 调 和 分 析 、 几 何 测 度 论 中 的 应 用 。Bouligand 1928 
年 引入 了 Bouligand 维 数 ,Poutrjagin 与 Schnirelman 在 1932 年 引入 了 覆盖 维 数 ， TE 
哥 洛 夫 (Kolmogorov) 与 Tikomirov F 1959 ££ SA T HZ 列 维 (Levy) 在 下 面 两 方 
面 的 研究 成 果 至 今 看 来 仍然 显得 非常 重要 : 其 一 ， 他 第 一 个 系统 地 研究 了 目 相 似 
集 ， 现 今 人 们 所 知 许多 自 相 似 集 的 性 质 几 乎 都 可 以 追溯 到 他 的 工作 ;其 二 ， 他 建立 
了 分 数 布朗 运动 的 理论 ， 可 以 说 ， 他 是 随机 分 形 理 论 研究 最 重要 的 先驱 者 之 一 。 以 
Salem 与 Kahane 为 代表 的 法 国学 派 则 从 稀薄 集 的 研究 出 发 ， 对 集 及 做 了 系统 的 研 
究 ， 体 现 出 了 其 独特 的 理论 方法 和 技巧 ， 并 在 调和 分 析 中 找到 了 其 重要 的 应 用 。 与 
此 同时 ， 维 数 的 乘积 理论 、 网 格 测度 计算 技巧 、 位 势 方法 及 投影 理论 等 均 先 后 被 建 
立 起 来 。 所 有 这 些 已 经 使 得 分 形 科 学 的 研究 具有 了 自己 的 分 析 特 色 和 方法 。 

尽管 这 一 阶段 的 研究 取得 了 这 么 多 成 果 ， 但 是 绝 大 部 分 人 主要 局 限于 对 数学 理 
论 的 研究 而 未 和 其 它 学 科 发 生 联系 。 另 一 方面 ， 物 理 、 地 质 、 工 程 学 等 学 科 则 已 产 
生 了 大 量 和 分 形 有 关 的 问题 ， 妃 切 需 要 新 的 理论 思想 和 有 力 的 工具 来 处 理 。 正 是 在 
这 样 的 背景 情形 下 ， 自 六 十 年 代 以 来 ， 曼 德尔 布 罗 特 以 其 独特 的 思想 ， 先 后 系统 深 


BENED a  : 


入 地 研究 了 海岸 线 的 结构 、 月 球 表面 、 地 貌 生 成 的 几何 性 质 、 强 噪音 干扰 下 的 电子 
通讯 等 自然 界 中 大 量 典型 的 分 形 现象 。 

第 三 阶段 て 197 年 至 今 ) ， 分 形 理论 在 各 个 领域 的 应 用 取得 了 全 面 的 发 展 ， 并 
形成 了 独立 的 学 科 。 曼 德 布 罗 特 将 前 人 的 结果 加 以 总 结 ， 集 其 大 成 ， 于 1975 年 以 
“分 形 ， 形 状 、 机 遇 和 维 数 ” 为 名 发 表 了 他 的 划时代 的 著作 。 他 是 第 一 个 跳出 传统 
物理 学 和 几何 学 的 人 。 在 这 本 专著 中 ， 他 第 一 次 系统 地 阐述 了 分 形 几 何 的 思想 、 内 
容 、 意 义 和 方 法 ， 标 志 着 分 形 几 何 作为 一 门 独立 的 学 科 自 式 诞生 。 但 是 ， 关 于 分 形 
是 否 是 一 门 数学 科学 的 讨论 在 这 个 时 候 也 此 起 彼 伏 ， 更 加 激烈 。 特 别 一 提 的 是 ， 
1988 年 以 来 B.B.Mandelbrot 5j S.Krantz 中 一直 在 为 分 形 的 价值 而 争论 不 休 。Krantz 
认为 “分 形 没有 明确 的 定义 ， 作 为 一 个 数学 家 ， 我 觉得 这 不 是 一 个 好 兆头 ”，“ 他 
们 产生 图 形 是 为 了 得 到 更 多 的 图 形 ， 而 不 是 为 了 得 到 更 深刻 的 思想 ”。Mandelbrot 
则 认为 “有 关 分 形 的 研究 加 深 了 我 们 对 自然 的 理解 。 如 果 我 只 是 证 明了 少数 几 个 定 
理 的 话 ， 那 么 用 这 些 定理 就 很 难 发 现 现在 还 没有 创立 的 或 潜在 的 研究 领域 ”。 这 场 
Mandelbrot 与 Krantz 最 具 代 表 性 的 争论 甚至 被 称 为 “是 一 场 关 于 数学 灵魂 的 战 
RU, 

争论 总 归 争 论 。 自 1975 年 以 来 ， 分 形 理论 在 数学 基础 和 应 用 方面 都 取得 了 相当 
迅速 的 发 展 。 在 物理 的 相 变 理论 、 材 料 的 结构 、 力 学 中 的 断裂 与 破坏 、 高 分 子 链 的 
聚合 、 酶 的 生长 机 理 研究 、 自 然 图 形 的 模拟 和 模式 识别 等 领域 取得 了 令 人 满意 的 成 
就 户 由 jos 。 近 十 几 年 来 ， 在 非 线性 应 用 学 科 和 计算 机 制图 的 刺激 和 推动 下 ， 分 形 
的 数学 理论 也 得 到 了 更 加 深入 的 发 展 ， 这 主要 体现 在 以 下 几 个 方面 ; 分 形 维 数 的 信 
计 与 算法 ， 分 形 集 的 生成 结构 ， 分 形 的 随机 理论 ， 动力 系统 的 吸引 子 理论 ， 分 形 的 
局 部 结构 ， 小 波 分 析 方 法 开始 在 分 形 分 析 中 得 到 应 用 。 这 其 中 特别 值得 一 提 的 是 ， 
BarnsieyPJeij 利 用 一 类 设 定 的 迭代 函数 系统 来 逼近 分 形 的 不 变性 ， 从 而 成 为 开创 了 
分 形 反 演 问题 的 先驱 。 后 来 ， 在 A.Arneodo、Z.R.Struzik、Vrscay、Forte、 朱 治 车 和 
李 后 强 等 人 PP 的 努力 下 ， 再 加 上 客观 上 图 形 压 缩 的 迫切 需要 ， 分 形 反 演 问 题 
的 研究 增 掀起 了 一 场 高 潮 。 在 此 期 间 ， 自 然 科学 中 的 分 形 学 术 论 文 旺 指数 增长 趋 
势 ， 哲学 社会 科学 领域 中 涉及 分 形 的 论文 和 书籍 也 不 断 增 加 。 国 内 国际 有 关 分 形 的 
专题 会 议 有 增 无 减 ， 特别 是 在 八 十 年 代 后 期 ， 令 人 感到 了 雷霆 万 钧 乙 势 。 国 际 学 术 
刊物 《混沌 、 孤 子 和 分 形 》(Chaos,Soljitons and Fractals,Pergamon Press,1991) 和 《分 
形 》 (Fractals,World Scientific,1993) 已 先后 问世 ， 分 形 学 科 从 此 就 有 了 属于 目 己 的 一 
块 阵地 。 
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“分 形 ” 一 词 传 入 我 国 以 后 ， 立即 在 各 行 各 业 中 引起 了 广泛 的 传播 ， 研 究 分 形 
理论 及 应 用 的 队伍 日 益 壮大 。 我 国政 府 和 学 术 界 很 早 时 候 就 都 对 分 形 理论 一 直 持 积 
极 的 肯定 态度 ， 并 给 予 了 极 大 的 支持 。 在 国家 863 攀登 计划 《 非 线性 科学 》 项 目 
中 ， 列 出 了 “分 形 的 数学 理论 ”《【 随 机 多 分 形 的 数学 理论 、 维 数 的 计算 方法 等 ) 、 
“分 形 的 物理 机 理 ” (多 重 分 形 结构 、 分 形 统计 模型 的 相 变 、 酶 结构 分 析 、 动 力学 
集团 的 生长 等 ) 两 大 发 展 方 向 ， 国 家 自然 科学 基金 申请 指南 中 也 已 列 出 “分 形 理论 
及 其 应 用 ”(A01020405) 内 容 。 在 分 形 的 数学 理论 发 展 方面 ， 国 内 南京 大 学 学 者 苏 维 
HAMA 1990 年 研究 了 分 形 与 局 部 紧 群 上 的 调和 分 析 ， 熊 金城 、 Ra. HALEY 
等 人 ?4 下 则 分 别 在 1994 年 至 1995 年 各 自 研 究 了 一 类 特殊 集合 的 维 数 问题 ， 张 
永 平 虽 在 分 形 微 积分 方面 提出 了 采用 多 项 式 最 佳 平方 逼近 来 定义 分 形 函 数 拟 导数 概 
念 的 方法 ， 胡 迪 稚 四则 在 研究 随机 分 形 方面 获得 了 较 大 的 成 功 。 朱 治 军 和 李 后 强 两 
人 国 对 分 形 反问 题 进行 了 较 深入 的 研究 ， 弥 补 了 国内 在 此 问题 研究 方面 上 的 不 足 。 


1. 1. 2 ”分形 及 其 维 数 的 定义 


什么 是 分 形 ? 事实 上 ， 到 目前 为 止 ， 分 形 还 没有 严格 的 数学 定义 ， 只 能 给 出 一 
般 描述 性 的 定义 。 有 关 分 形 的 概念 ， 科 学 史上 有 记载 的 最 早 是 由 Hausdorff 于 1919 
年 引入 ,随后 经 Besicovitch 于 1935 年 和 Mandelbrot 于 1975 年 加 以 改进 和 发 展 。 
Mandelbrot 是 美 籍 法 国 科 学 家 ， 在 撰写 他 的 分 形 学 莫 基 性 著作 《分 形 ， 机 遇 和 
维 数 》 中 碰 到 法 语 动词 frangere (EI) 变 来 的 形容 词 fractus， 联 想到 英语 中 的 同 源 
ial fracture (断裂 ) 和 形容 词 fraction (分 数 ) , 他 在 这 些 词 的 基础 上 创造 了 分 形 
fractal 一 词 。Fractal 最 初 的 意思 是 破碎 的 和 不 规则 的 ， 并 将 分 形 定义 为 具有 某 中 意义 
下 的 自 相 似 集合 。Mandelbrot 也 曾 对 分 形 做 了 一 个 尝试 性 的 定量 刻画 ， 认 为 分 形 是 
Hausdorff 维 数 严格 大 于 其 拓扑 维 数 的 集合 ， 但 这 些 定义 都 不 够 全 面 和 精确 。 英 国 数 
学 家 在 其 著作 中 认为 ， 分 形 的 定义 应 该 以 生物 学 家 定义 生命 的 类 似 方式 给 出 ， 也 碘 
是 说 不 寻求 分 形 的 严格 确切 简明 的 定义 ， 而 是 寻求 分 形 的 特性 ， 在 这 种 观点 的 指导 
下 ， 人 们 将 分 形 看 成 是 有 具有 如 下 性 质 的 一 类 集合 : 
1. 具有 精细 的 结构 ， 即 在 任意 小 的 比例 尺度 内 包含 了 整体 ， 类 似 于 生物 中 的 全 
Rf 
2. 分 形 既 可 以 是 几何 图 形 ， 也 可 以 是 由 “功能 ”或 “信息 ”建立 起 来 的 一 种 数 
理 模 型 ， 
3. 共有 不 规则 性 以 致 它 的 整体 与 与 局 部 都 不 能 用 传统 的 几何 语言 来 表达 ; 
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4， 具 有 某 种 自 相似 的 形式 ， 可 以 是 严格 自 相 似 的 也 可 以 是 统计 意义 上 的 相似 
性 ， 自 然 界 中 的 大 部 分 分 形 均 属于 后 一 类 ; 

5， 可 以 同时 具有 形态 、 功 能 和 信息 三 方面 的 自 相似 性 ， 也 可 以 只 是 其 中 某 一 方 
面 的 自 相似 形 

6 相似 形 具 有 层次 结构 方面 的 差异 。 数 学 理论 上 的 分 形 具 有 无 限 符 套 的 层 钦 结 
构 ， 而 自然 界 中 的 分 形 只 是 有 有 限 层 次 的 媒 套 ， 在 进入 到 一 定 的 层次 以 后 才 
会 有 分 形 的 規律 

7. 相似 形 具 有 级 别 上 的 差异 〈 生 成 元 的 使 用 次 数 或 者 放大 倍数 ) 。 AIR 
近 ， 则 愈 相 似 。 级 别 相差 愈 大 ， 相 似 形 也 就 愈 差 。 整 体 的 级 别 最 高 ， 最 低 的 
是 级 生成 元 。 可 以 用 无 标 度 区 间或 者 标 度 不 变形 范围 加 以 表示 。 目 然 界 的 分 
形 往往 均 具 有 一 个 最 大 和 最 小 的 标 度 ， 在 无 标 度 区 域内 ， 对 象 才 具有 分 形 规 
律 ， 否 则 ， 一 旦 越过 无 标 度 区 ， 自 相似 性 也 就 随 之 消失 ; 

8. 一 般 情况 下 ， 以 某 种 方式 定义 的 分 形 维 数 大 于 该 分 形 的 拓扑 维 数 ; 

9. 在 多 数 令 人 感 兴趣 的 情况 下 ， 分 形 常常 以 非常 简单 的 方法 来 定义 ， 或 许 是 一 
种 递归 方式 ; 

国内 学 者 也 对 分 形 的 定义 进行 了 专门 的 研究 ， 吴 敏 金 乓 在 其 分 形 信息 导论 一 书 
中 讨论 了 分 形 的 定义 问题 ， 认为 分 形 是 一 个 不 断 分 化 演变 的 动态 过 程 ， 并 从 演化 的 
观点 在 这 方面 做 出 了 许多 开创 性 的 工作 。 

分 形 度 量 则 是 另外 一 个 与 分 形 概念 相关 的 比较 敏感 的 问题 ， 大 多数 分 形 集 的 定 
量 刻 画 和 判别 要 借助 于 参数 一 一 维 数 概念 ， 因此 维 数 的 定义 也 是 一 个 人 们 比较 关心 
的 问题 。 

历史 上 对 维 数 的 认真 研究 似乎 直到 19 世纪 后 半 叶 才 开 始 。 此 前 ， 虽 然 Peano 和 
Minkowski 都 对 维 数 的 定义 提出 过 重要 的 思想 ， 但 都 只 是 直观 上 的 描述 ， 还 缺乏 严格 
的 数学 定义 。 最 初 人 们 党 把 维 数理 解 为 所 研究 几何 对 象 的 坐标 数 ， 是 描写 空间 内 一 
个 点 所 需 的 最 少 实 参数 的 个 数 。 在 物理 世界 里 ， 自然 界 客体 的 维 数 与 所 在 空间 的 维 
数 有 时 是 不 同 的 客体 的 维 数 与 测量 尺度 有 关 。 从 远 处 观察 一 个 毛线 团 ， 人 们 可 把 它 
看 成 是 一 个 点 ， 而 从 近 处 观察 它 ， 则 人 们 会 认为 它 是 一 条 空间 曲线 。 经 典 的 Peano 
曲线 例子 就 是 通过 连续 映射 将 一 个 区 间 映 射 到 一 正方 形 ， 在 这 样 的 情形 下 ， ERE 
数 观 点 也 就 失去 了 支持 。 维 数 的 一 个 重要 进展 应 该 首 推 Hausdorff 于 1919 年 引入 的 
维 数 概 念 。 该 维 数 以 Hausdor 人 f 测 度 为 基础 ， 且 适合 于 任何 集合 ， 从 定义 出 发 作 维 数 
的 计算 是 相当 困难 的 ， 这 种 计算 上 的 困难 也 极 大 地 限制 了 Hausdorff 维 数 的 应 用 。 但 
要 了 解 分 形 , 人 们 就 有 必要 了 解 Hausdorff 测度 与 Hausdorff 维 数 。 随 着 分 形 应 用 研 
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究 的 不断 深 入 , 分 形 维 数 的 定义 多 种 多 样 ， 只 能 根据 不 同 的 研究 对 象 来 使 用 各 别 的 
计算 方法 。 除 Hausdorff 维 数 以 外 ,还 有 相似 维 数 、 计 盒 维 数 、 容 量 维 数 、 填 充 维 数 、 
Kolmogorov 维 数 、 Lyapunov 第 数 等 Gi Ze [8] ALANS 6-15 $8 "1010415 , 而且 这 些 维 
数 之 间 的 关系 还 没有 理 顺 。 大 多 数 分 形 维 数 的 定义 均 是 基于 “尺度 5 下 的 度量 ”这 
一 思想 。 设 下 是 一 分 形 集 ， 对 于 每 个 6 > 0 ， 忽 略 尺 度 小 于 的 不 规则 性 ， 并 且 考 虑 
MEHM, (FES 2 0 时 的 情形 。 若 存在 两 个 常数 k>0 和 s >0， 使 得 MM(F) 满 足下 
E BB: 
M ,(F)-kó^ 

则 称 集合 Ff 具有 维 数 s， 而 可 作为 集合 在 s 维 下 的 长 度 。 

另外 一 种 关于 分 形 维 数 的 理解 是 由 我 国学 者 赵 辉 和 孙 博 文 吧 提出 来 的 。 他 们 认 
为 维 数 是 集合 层次 的 量 值 标号 。 当 人 们 从 不 同 的 角度 考 得 一 个 集合 中 的 元 素 时 ， 残 
会 发 现 这 些 元 素 具 有 差异 性 。 此 时 ， 我 们 可 以 赋予 集合 某 种 “层次 结构 ”， 按 这 个 
层次 结构 对 集合 的 元 素 去 分 类 处 理 ， 每 一 类 中 的 元 素 属于 同一 层次 ， 不 同 层 次 中 的 
元 素 相 对 于 考查 内 容 来 说 具有 质 的 不 同 。 在 确定 的 层次 中 去 考察 不 属于 该 层次 的 元 
素 时 ，“ 奇 异性 ”也 就 难免 产生 了 。 例 如 ， 在 一 维 层次 中 【使 用 测度 为 1〉 去 观测 一 
条 有 界 直 线 的 测度 时 ， 得 到 的 是 一 个 有 界 实数 ， 同样 去 观测 Von.Koch 曲线 的 测度 
时 ， 则 得 到 数学 上 的 “无 穷 大 ”结果 。 相 反 ， 在 Von.Koch 曲线 层次 中 《使 用 测度 为 
d 观测 一 条 直线 时 ， 得 到 的 是 “无 穷 小 ”的 结果 ， 而 去 观测 Von.Koch Hee 
则 都 是 一 个 有 界 实 数 。 由 此 可 见 ， 集合 的 奇异 性 或 复杂 性 的 产生 是 由 于 人 们 跨越 层 
次 考察 事物 导致 的 。 

综 上 上 所 述 ， 由 于 分 形 的 确切 定义 疝 未 给 出 ， 因 此 分 形 的 度量 一 一 分 维 数 也 就 难 
以 精确 定义 ， 人 至 今 仍 缺乏 统一 的 计算 公式 与 算法 。 


1. 1. 3 目前 分 形 几何 学 的 研究 工具 和 研究 方法 


分 形 集 无 论 在 局 部 还 是 在 整体 都 是 不 光滑 的 ， 这 就 决定 了 经 典 数 学 中 与 微 积分 
有 关 的 许多 分 析 工 具 在 这 里 不 能 发 挥 其 作用 。 分 形 早 期 的 主要 研究 工具 是 点 集 拓扑 
与 集合 测度 论 ， 这 主要 体现 在 Besicovitch 和 Hausdorff HAMMAR LE. Ab 
种 重要 的 工具 是 Julia 等 人 所 采用 的 复 动力 系统 理论 。 借 助 计算 机 的 强大 能 力 ，Julia 
集 和 Mandelbrot 集 鸭 的 精细 复杂 而 美丽 诱 人 的 结构 展示 于 世人 的 面前 ， 这 也 是 分 形 
能 够 广泛 而 迅速 为 人 们 所 接受 的 一 个 重要 原因 。 


ee 


白 从 八 十 年 代 后 期 MF.BarnsleyDJ 四 采用 选 代 函 数 系统 IFS(Iterated Funtion 
System) 构 造 出 许多 分 形 以 来 ， 迭代 函数 理论 的 研究 方法 和 手段 也 就 组 成 了 研究 分 形 
的 重要 部 分 , REPARER FIERE. 由 于 单位 区 间 上 每 一 点 都 可 以 看 
做 是 一 个 十 进 制 的 符号 序列 ， 因 此 ， 对 于 和 一 般 符号 序列 复杂 性 研究 有 关 的 连 分 数 
理论 、 遍 历 理 论 、 符 号 动力 系统 理论 都 是 现 阶段 研究 分 形 的 新 工具 。 

特别 值得 一 提 的 是 ， 近 几 年 来 ， 有 “数学 显微镜 ” 之 称 的 小 波 分 析 理 论 为 具有 
“无 穷 精细 结构 ”的 分 形 的 研究 提供 了 一 个 绝 好 的 分 析 工 上 基 ， 通过 小 波 变换 ， 人 们 
可 以 看 到 分 形 对 象 的 丰富 细节 ， 为 推测 分 形 生成 的 动力 学 根源 提供 了 方便 。 这 主要 
体现 在 Argoul、Areodo、Holschneider、Struzik、 李 后 age A Pub 的 工作 中 。 另 
外 ， 本 人 在 论文 中 也 做 了 一 些 有 关 这 方面 的 研究 工作 。 

在 方法 上 ， 除 理论 性 的 研究 外 ， 通 过 计算 机 强大 的 计算 能 力 和 图 形 能 力 ， 开展 
“实验 数学 ”方法 在 分 形 的 研究 中 也 将 起 到 重要 的 作用 。 在 实验 的 基础 上 进行 理论 
的 思考 ， 实 验 有 助 于 人 们 更 好 的 借助 “直觉 ”理解 数学 规律 ，Mandelbrot 在 倒立 分 
形 几何 时 就 曾经 使 用 该 方法 。 美 国 大 学 已 经 设立 了 实验 数学 系 。 随 着 计算 机 辅助 方 
法 的 大 量 使 用 和 推广 ， 有 理由 相信 ， 实验 数学 有 可 能 成 为 一 门 探索 性 的 独立 学 科 。 


1. 1. 4 目前 分 形 理论 研究 中 存在 的 主要 问题 及 动 问 


分 形 理论 是 在 发 展 中 不 断 完善 和 走向 成 熟 的 。1993 年 在 中 国 科 技 大 学 召开 的 全 
国 第 三 界 分 形 学 及 其 应 用 研讨 会 上 ， 李 后 强 等 人 Pei 就 此 后 一 个 时 期 分 形 学 及 其 应 
用 研究 提出 了 以 下 5 个 大 问题 和 10 个 方面 。 这 些 问 题 在 一 定 程度 上 代表 了 目前 分 形 


研究 的 热点 。 


1， 如 何 判断 一 个 对 象 是 分 形 或 多 分 形 的 问题 

如 前 所 说 ， 分 形 没 有 完整 而 精确 、 一 致 公认 的 定义 ， 要 判断 分 形 与 多 分 形 是 很 
困难 的 。 也 有 人 认为 ， 分 形 就 象 生命 一 样 理 解 其 含义 就 可 以 了 。 在 理论 上 ， 上 日 仿 身 
性 是 分 形 产生 的 重要 方式 ， 所 以 有 关 自 仿 射 人 性 方面 的 研究 应 该 是 一 个 重点 方面 。 但 
对 于 实际 中 出 现 的 分 形 进行 判断 依然 是 一 个 复杂 而 未 解 的 问题 。 


2. 分 数 维 所 代表 的 实际 意义 问题 

分 数 维 到 底 有 什么 用 ? 这 个 问题 是 关乎 分 形 这 门 学 科 的 优势 及 其 存在 价值 的 六 
事 。 在 这 方面 Hausdorff 维 数 的 意义 比较 明确 ， 它 可 定量 地 表示 一 个 集合 规则 或 不 规 
则 的 的 集合 尺度 ， 而 在 动力 系统 方面 ， 则 大 体 上 表示 独立 变量 的 个 数 。 多 分 形 的 奇 
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异 谱 主要 表征 多 分 形 的 奇异 性 。 在 时 间 序 列 中 所 用 的 关联 维 数 或 广义 维 数 也 反映 了 
序列 的 某 种 相关 性 和 复杂 程度 。 继 续 探索 各 类 维 数 的 实际 意义 仍然 是 今后 工作 的 一 
^P. 


3. 分 形 的 动力 学 机 制 问题 

研究 分 形体 集团 生长 的 时 间 演 化 规律 和 结构 的 标 度 行为 ， 可 以 改变 分 形 理论 中 
仅仅 只 对 分 形 进行 形态 措 述 的 单调 而 被 动 的 局 面 。 非 平衡 物理 学 中 的 耗 散 结 构 理 论 
和 协同 学 中 许多 概念 有 必要 引入 分 形 的 研究 领域 ， 另 外 对 于 分 数 阶 微分 方程 、 非 线 
性 发 展 方程 、 神 经 网 络 方面 也 应 该 予以 一 定 程度 的 关注 。 要 重视 随机 噪音 对 分 形 系 
统 的 影响 。 这 方面 最 漂亮 的 成 果 是 ， 在 对 扩散 凝聚 集团 DLA 生长 过 程 的 研究 观测 
中 , 国外 Witten 和 Sander 成 功 地 建立 了 其 分 形 生 长 动力 学 模型 。 


4. 分 形 的 重 构 问 题 

广义 而 言 ， 分 形 重 构 是 指 一 个 被 认为 是 分 形 的 图 形 能 否 以 某 种 指定 的 方式 生成 
它 。 这 是 动力 系统 研究 领域 中 的 反问 题 。 分 形 所 表现 出 的 自 相 似 性 和 拟 自 相似 性 基 
本 上 都 和 迭代 、 递 归 过 程 有 关联 。Barnsley 和 Demkol 中 最 早 进行 了 用 不 变 测度 IFS 
逼近 目标 测度 的 工作 。 抢 量 法 和 拼 贴 定理 Eu] 是 研究 该 类 问题 的 两 个 基本 方法 。 这 
主要 体现 在 Vrscay 和 Forte 等 人 的 一 系列 研究 工作 中 。 近 几 年 来 ， 把 小 波 分 析 方 法 
应 用 于 分 形 反 问题 研究 ， 也 获得 了 许多 好 的 结果 《可 参阅 本 论文 第 三 章 内 容 ) 。 分 
形 重 构 问 题 对 于 图 形 和 数据 压缩 有 着 非常 重要 的 意义 ， 这 方面 研究 工作 急需 继续 努 
力 和 深化 。 


5. 关于 Julia 集 和 Mandelbrot 集 的 问题 

动力 系统 中 的 Julia 集 和 Mandelbrot 集 是 复 多 项 式 迭 代 后 的 结果 。 研 究 各 类 多 项 
式 〈 如 整 系数 多 项 式 和 代数 整 系数 多 项 式 等 ) 的 填充 Julia 集 内 准 周期 点 的 代数 特 
性 ， 有 关 Mandelbrot 集 是 否 是 局 部 连通 的 及 边界 维 数 的 计算 ,这 些 问 题 一 直 是 分 形 研 
究 中 的 难点 和 热点 之 一 。 


以下 10 个 方面 有 待 尽 快 研究 以 取得 突破 : 
1. 随机 多 分 形 的 数学 理论 ; 
2. 分 形 曲 线 的 导数 问题 ; 
分 形 集 往往 存在 分 数 阶 导数 ， 这 种 分 数 阶 导数 对 刻画 分 形 对 象 有 重要 的 作 


の 
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用 并 由 此 建立 一 个 完整 的 分 数 阶 微 积 分 体系 。 
.分 形 的 小 波 分 析 及 小 波 变 换 产生 分 形 的 关系 问题 ; 
.分 形 分 维 数 的 计算 问题 ; 

.多 分 形 的 热力 学 及 相 变 普 适 性 判 据 问 题 
原子 、 分 子 的 分 形 问 题 : 

. 胖 分 形 及 重 正 化 泥 沌 问题 ; 

. Bicis ESL AR TURPE IRL; 

， 生 物 膜 的 分 形 结构 及 其 与 膜 病 变 的 关系 问题 ; 
10， 图 形 的 分 形 压缩 问题 ; 


On A UC 


L2 ”小 波 理论 及 其 在 分 形 有 关 问 题 中 的 应 用 
1. 2. 1 小 波 的 发 展 及 定义 


小 波 变换 分 析 方 法 是 八 十 年 代 中 后 期 逐渐 发 展 起 来 的 一 种 数学 分 析 方 法 ， 是 调 
和 分 析 近 四 十 年 来 的 一 个 突破 性 进展 。 它 在 保留 了 Fourier 分 析 优 点 的 同时 ， 又 弥补 
了 Fourier 分 析 不 能 作 局 部 及 非 平稳 分 析 等 不 足 ， 在 理论 分 析 和 实际 应 用 中 得 到 乏 亏 
的 发 展 ， 目 前 正 受 到 数学 界 和 工程 界 的 极 大 重视 SIP] 。 

小 波 分 析 是 数学 家 、 应 用 数学 家 和 信号 处 理 、 数 据 处 理工 程 师 们 在 各 生 工 作 领 
域 分 别 独立 发 现 的 。 小 波 分 析 方 法 的 提出 ， 最 早 应 追溯 于 Haar 1910 年 提出 的 规 
范 正 交 基 。1938 年 ,Littlewood-Paley 对 Fourier 级 数 建立 了 L-P 理论 ， 即 按 二 进 制 频 
率 成 分 分 组 ，Fourier 变换 的 相位 变换 本 质 上 不 会 影响 函数 的 形状 与 大 小 。 在 1965 年 
Galderon 发 现 了 再 生 公 式 , 它 的 离散 形式 已 接近 小 波 展开 ， 只 是 还 无 法 得 到 一 组 正 交 
RWE. Stromberg?! TE 1981 年 对 Haar 系 进行 了 改进 ,证 明了 小 波 函 数 的 存在 性 。 
1982 年 Battle 在 构造 量子 场 论 中 采用 了 类 似 于 Galderon 再 生 公 式 的 展开 式 。 

小 波 概念 的 真正 出 现 应 该 算 于 1984 年。 和 Fourier 分 析 一 样 ， 小 波 分 析 包 括 两 
个 部 分 的 内 容 ， 小 波 变 换 和 小 波 级 数 。 

小 波 变 换 是 由 法 国 科学 家 J.Morlet1984 年 的 在 分 析 地 震波 的 局 部 特性 时 引入 
的 ， 并 首次 提出 小 波 (Wavelet) 这 一 术语 。 所 谓 小 波 就 是 小 的 波形 ，“ 小 ” 指 其 具有 
衰减 性 ，“ 波 ” 指 其 波动 性 ， 具 有 振幅 正 负 相 间 的 震荡 形式 。 随 后 ，J.Morlet 和 物理 
学 家 Grossman 等 人 共同 建立 了 完整 的 连续 小 波 变换 几何 体系 ， 其 基础 是 平移 和 伸缩 
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下 的 不 变性 ， 这 使 得 能 将 一 个 信号 分 解 成 对 空间 和 时 间 尺 度 的 独立 贡献 ， 同 时 义 个 
丢失 原 有 信和 号 的 信息 。 

小 波 级 数理 论 则 依赖 于 小 波 基 的 发 展 。 法 国 马 赛 大 学 的 数学 家 Y.Meyer VH 在 
1985 年 将 J.Morlet 的 思想 与 数学 研究 中 已 存在 的 奇异 积分 算 子 理论 联系 起 来 ， 构 井 
了 具有 一 定 衰减 性 的 光滑 函数 ， 其 二 进 制 尺度 位 移 体 系 构成 了 忆 (R) 空 间 的 一 组 标准 
正 交 基 , 同時 Y.Meyer 与 A.Grossman 及 比利时 女 数 学 家 I.Daubechies Biss] 4 He, 通 
过 构造 2(R) 上 的 一 准 正 交 完备 集 的 方法 选取 了 小 波 空间 的 离散 子 集 一 一 框架 ， 并 证 
明了 一 维 小 波 函 数 的 存在 性 。 稍 后 ，P.G.Lemarie 则 将 这 一 结果 推广 到 n 维 情况 ， 并 
和 G.Battle 给 出 了 具有 指数 衰减 的 小 波 函 数 。I.Daubechies 4711988 年 发 表 了 她 的 著 
名 长 篇 论文 ， 证 明了 具有 有 限 紧 子 集 的 正 交 小 波 基 的 存在 性 。S.Mallat Pre 在 
1989 年 将 计算 机 视觉 领域 内 的 多 尺度 分 析 思 想 引 入 到 小 波 分 析 中 ， 提 出 多 分 辨 率 分 
析 MRA(Multiresolution Analysis) 概 念 ， 统 一 了 在 此 之 前 所 有 具体 正 交 小 波 的 构造 ， 
并 提出 了 相应 的 分 解 和 重 构 快 速算 法 ， 其 地 位 相当 于 FFT 在 经 典 Fourier 分 析 中 的 地 
位 。1988 年 ，LDaubechies 在 美国 NSF/CBMS 主办 的 小 波 专 题 研讨 会 上 进行 了 10 次 
讲演 四， 引起 了 广大 数学 家 、 观 察 学 家 、 物 理学 家 甚至 某 些 企业 家 的 高 度 重视 ， 由 
此 将 小 波 分 析 的 理论 发 展 与 实际 应 用 推 向 了 一 个 高 潮 。 

通常 小 波 有 如 下 一 些 定义 : 
定义 1.2.1 (Grossman and Morlet) 如 果 函 数 y(x)e L (R)I Fourier Æ% p(o) WE 
下 式 要求 


则 称 y(x) 为 三 (RR) 上 的 一 个 小 波 函 数 。 

定义 1.2.2 (Littlewood-Paley-Stein) 如 果 函 数 y(x) e ど (人 W) 的 Fourier 変換 (の ) 満 足 
下 式 

y vp" o) = l(a.e) 

別称 w(x) e P (R) AMERS. 

定义 1.2.3 (Stromberg and Franklin) P? 如果 na v(2'x-k) Uke Z)o P (m)iti— 
组 规范 正 交 基 ， 则 称 y(x)e D(9)7s— ARR. 

定义 1.2.4(Y.Meyer) 如 果 一 个 实 变量 函数 ywy(x) 满 足以 下 四 个 性 质 ， 则 称 y(x) 为 一 
个 m(m eZ) 美 的 小波 : 

(1) Hm=0, AA v(x)er(9*): 若 m>1， 则 w(x) 及 它 的 直到 m 阶 导数 都 属 

于 (8)。 
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(2) y(x) REM HEI m E SMAI LEME (局部 正則 性 ) 。 
(3) 对 于 0<k<m， 有 [xyw(z) な =0 (消失 知 条件 ) 。 


(4) 2” y(2’x—k) (j,ke Z) dier P (9) S d iE. 
小 波 函 数 的 存在 证 明 〔 构 造 方法 ) 以 及 各 种 小 波 算法 都 和 多 分 凑 分 析 MRA AUS 
度 函 数 有 着 密切 的 关系 。 为 此 下 面 引 入 多 分 辨 分 析 和 尺度 函数 的 概念 。 
设 亿 } ELON TEP, BY} BERN TEAM, 如果 
满足 如 下 条 件 : 
(1) V, cV. NV, ={0}. Ur, =P R) 
(2) fer, c f(2x)eV,,, 
(3) f(x)e y, f(x-k)e V. keZ 
(4) 存在 函数 g(x)e Vo» 使 得 {g(x - fez FJ KV, 的 Riesz 基 
in RRR (x) e 了 (RR) 能 产生 二 (8%) 的 一 个 多 分 辨 分 析 ， 而 EH? ,所 张 成 的 
线性 空间 ， 即 就 是 
V, = clos d, € Z) 
b,, =27 (o xk) ke z 
则 称 g(x) 为 三 (RR) 的 一 个 尺度 函数 。 然后 ， 借 助 尺度 函数 g(x) 就 可 以 构造 出 小 波 函 
数 y(x) 来 。 
存在 尺度 系 数 je//, plaire iar 
ie) = Dhglor 8), bez er 


4 v(x)- Cy h, 92x - k) 


W, = poly (5) - 2^ 2 volt) ez} 
则 有 W, LV, W, OV, =V; AML = @W,» WL 
Wah pep 5L Qm 中 的 一 个 标准 正 交 小 波 基 。 
対 手 母 小波 函数 y(x), 当 x 、a、5 连续 变换 时 ， 进 行 伸缩 和 平移 后 形成 的 一 族 
小 波 记 为 fy,,}， 则 有 连续 小 波 
a) = Avi M). aben, a>0 
对 任何 f(x)e 2 (9), n w 人 xz] 作 为 窗 函数 的 连续 小 波 变 换 定义 为 
wat) =z [rev —* la 
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可 见 小 波 变 换 已 将 f(x) BO, 中 心 在 x= ヵ , REA 2A, BE) f(x) 的 時 一 
频 局 部 化 。x 一 w 平 面 上 的 时 频 窗口 为 
[5— ~ady, bean] - Lag Tag 


通常 取 @。 > 0， 由 上 式 可 见 ， 在 高 频 时 候 ， 时 窗 变 帘 ;， 在 低频 时 候 ， 时 窗 变 
宽 。 这 就 是 小 波 变换 所 特有 变焦 距 性 质 。 这 一 特性 克服 了 窗口 傅立叶 变换 的 缺陷 ， 
从 而 使 它 在 突变 信号 处 理 上 具有 了 特殊 的 地 位 和 作用 。 

若 w(x) 満 足下 式 


(o) 
C,=| ———d 
cL e| 
则 称 为 y(x) 允许 小 波 。 此 时 ， v 
f(x)= | [LONE dadb 


-00-0 


oc 
切 相关 的 另 一 个 概念 一 一 框架 ， 它 是 规范 正 交 基 的 推广 "2 。 
定义 1.2.4 jeJ ( 指 本 集 ) H, ¢ eH (Hilbert 空间 ) 。 若 存在 两 个 正 数 A 
和 B, 使 得 下 式 成立 
AT sre] snl vre n 


Ws fo; e J) 为 互 的 一 个 框架 ，A,B 为 框架 的 界 。 若 A- B. RUE, e] 
一 个 紧 框 架 (Tight Frame). 


1. 2. 2 分 形 研究 中 常用 的 几 个 小 波 


1. Daubechies 紧 支撑 正 交 小 波 ( 基 )) hh 
具有 紧 支 集 的 正 交 小 波 基 在 数字 信号 的 小 波 分 解 过 程 中 可 以 提供 有 限 长 度 的 更 
实际 、 更 具体 的 数字 滤波 器 ， 而 且 ， 它 和 Haar 小波 (MAM BAH) 等 相 比 ， 
Daubechies 连续 的 紧 支 集 正 交 小波 可 以 有 任意 高 阶 消失 和 矩 。 对 于 工 阶 消失 矩 的 小 波 ， 
其 滤波 器 有 的 长 度 N = 25 ， 并 且 工 -1 次 连续 可 导 。 
由 MRA 理论 可 知 ， 尺 度 肖 数 和 小 波 函 数 均 满足 双 尺 度 方程 ; 
(x) = = V2 22 ん @(2x - n) 


ve) 4 ay (-1) 4... 0(2x - n) 
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为 了 保证 紧 支 性 ， 就 必须 上 式 右边 含有 有 限 项 ， 作 适当 的 平移 变换 ， 上 面 两 式 
就 可 写成 
(x)= 25 hax - n) 
v(x) - V2 Y gd n) 


nzl-N 


构造 上 式 中 紧 支 p(x) (等 价 与 寻找 特征 多 项 式 H(z) ae ES 


归结 如 下 儿 个 步骤 : 
(D 逃 定 一 正 整数 ん >2. 
(2) 选 定 一 多 项 式 R ， 使 它 满足 下 面 三 个 表达 式 


(ke 


P (y)= M NA 


j 
(3) 寻找 一 实 系 


=P, jr 
(4) 则 得 H(z)= k 224] Q(z). 


2. Morlet 小波 
Morlet 等 人 利用 Gauss 函数 进行 移 位 展 缩 ， 定 义 了 如 下 的 Morlet 小 波 函 数 
V, (w)= DU ー AA , VKeR 
足 ゆ 。(0)= 0, X Fourier 逆 变 换 就 是 满足 容许 条 件 的 小 波 
ve) eet Ve n 


在 上 迷 小 波 解 析 式 中 , AKAROA, 第 二 項 可以 省略 揮 。 
3. Mexican 帽子 小波 
帽子 小 波 是 Gauss 函数 的 第 二 阶 导 数 ， 不 仅 具 有 无 穷 阶 光滑 性 ， 而 且 具 有 指数 
衰减 率 。 
= の ojze", (0)=0 
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这 种 满足 容许 条 件 的 解析 小 波 具 有 时 域 和 频 域 良好 的 局 部 性 质 ， 常 被 用 来 研究 
分 形 的 局 部 性 质 。 


1. 2. 3 小 波 在 分 形 研究 中 的 应 用 


小 波 理论 中 的 多 分 辨 分 析 思 想 ， 体 现 了 人 们 认识 和 识别 形体 的 过 程 遵循 了 一 种 
从 低 分 辨 率 到 高 分 辨 率 的 原理 。 对 于 分 形 的 研究 观察 也 正 是 这 样 的 方式 。 从 大 到 小 
的 不 同 尺度 变换 ， 在 越 来 越 小 的 尺码 上 观察 越 来 越 丰富 的 分 形 细节 问题 。 小 波 分 析 
和 分 形 具 有 深刻 的 内 在 联系 ， 两 者 在 尺度 变换 上 具有 一 致 性 ,分 形 是 一 种 几何 语言 ， 
而 小 波 是 一 种 分 析 局 域 奇 异性 的 工具 ， 两 者 相得益彰 。 

与 此 同时 ， 小 波 在 信号 处 理 、 图 象 边缘 检测 等 的 1 领域 中 的 广泛 应 用 为 小 波 
在 分 形 领 域 的 应 用 提供 了 方法 和 经 验方 面 的 借鉴 。 


小 波 在 以 往 分 形 研究 中 的 应 用 主要 在 以 下 几 个 方面 进行 。 


l. 在 分 形 反 问题 方面 的 应 用 : 

在 分 形 反 问题 研究 中 ，1993 年 A.Arneodo 首先 利用 小 波 对 一 维 饼 格 (Cookie- 
Cutters) 进 行 了 分 析 获 得 成 功 ， 从 而 开创 了 小 波 在 这 一 方向 应 用 的 先河 ， 葛 定 了 分 形 
反 演 的 理论 基础 。 随 后 A.Ameodo 等 四 人 的 方法 又 在 DLA 角 向 集 的 研究 中 发 现 了 存 
在 着 Fibonacci 结构 。 小 波 处 理 反 问题 的 关键 是 特征 点 的 选择 问题 。 后 来 的 人 们 几乎 
均 是 围绕 着 这 一 思想 展开 的 。Arneodo 最 初 选择 小 波 变 换 初 始 极 大 值 点 作为 特征 后， 
而 Z.R.Struzik 的 则 选择 小 波 极 值 分 岔 点 为 特征 点 ， 我 国学 者 朱 治 军 和 李 后 强 P3 选 
取 最 末 极 大 值 点 作为 特征 点 ， 其 计算 量 要 较 前 两 种 少 ， 并 可 在 较 大 的 尺度 上 进行 判 
断 。 


2. 在 各 种 类 型 分 形 结构 分 析 中 的 应 用 : 
分 形 的 典型 性 质 是 在 小 标 度 下 表现 出 某 种 自 相 似 性 。 局 部 目 相似 性 也 就 意味 着 
分 形 测度 u EA x 附近 的 标 度 变 化 满足 下 式 : 
u(B(x,,A«)) ~ a* pl B(x, E) 
其 中 B(x, e) PE x, ACME, ox, ) 是 局部 酸度 指数 ( 或 局部 Holder 指 
数 ) 。 为 了 研究 其 局 部 奇异 性 ， 因 此 小 波 可 被 用 来 分 析 分 形 的 测度 。 


当 A ORT, 其 标 度 特征 为 : 
W, (Aa, x, + Ab) ~ Ano (a, x, +b) 
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这 表明 了 分 形 测度 4 具有 相同 指数 w(xo ) 

这 方面 的 工作 主要 包括 有 均匀 和 非 均 邹 Cantor 集 的 小 波 变 换 、 A.Arneodo 等 人 
对 受 限 扩散 凝聚 模型 DLA(The Diffusion-Limited Aggregation) 的 小 波 变换 、 单 标 度 和 
双 标 度 雪 花 模型 的 小 波 变换 等 。 

3. 在 随机 分 形 分 数 布朗 运动 或 者 区 过 程 》 方面 的 小 波 分 析 研究 : 

分 数 Brown 运动 FBM(Fractional Brownian Motion) 是 一 个 非 平稳 的 具有 和 零 均值 的 
Gauss 随机 函数 1” 池 3 器 ， 在 描述 自然 界 和 社会 经 济 生 活 中 发 生 的 许多 事件 现象 〈 诸 
如 电子 元 器 件 中 的 噪音 、Internet 网 络 流量 、 股 票 价格 的 涨 落 等 ) DUS Re 
要 的 作用 。 其 定义 如 下 : 

B, (0)=0 


B80) = ke-80 - Cs)" pa) (6-57 an) 


T(H 4 0.5) 
其 中 B(f) 是 通常 的 Brown 运动 ，0 <H «1. 4H =~ Rt, FBM 就 是 Brown 运动 。 


主要 是 利用 分 数 布朗 运动 函数 B, (7) 的 小 波 变 换 结果 分 析 该 随机 过 程 具有 的 统计 
自 相似 特 征 。 最 早 利用 小波 方 法 的 是 Minnesota 大 学 的 A.H.Tewfik"?), 他 1991 年 发 
现在 不 同 变换 尺度 上 离散 小 波 系数 具有 相关 性 ， 系数 自 相 关 和 互相 关 函 数 按 抛物 线 
以 比 自身 自 相关 还 要 快 的 数 率 衰减 ， 衰减 率 主要 受 控 于 分 析 小 波 消失 和 矩 的 大 小 。 几 
平 在 同一 年 , 法 国 Lyon 物理 研究 所 的 Patrick Flandrin toj 研究 了 FBM 模型 的 小 波 合 
成 问题 。 国 内 华中 理工 大 学 的 王晓军 等 人 何在 离散 FBM 场 理论 的 基础 上 阐述 了 分 形 
小 波多 分 辨 图 像 识 别 的 概念 ， 提 出 了 一 种 高 效 快速 的 识别 算法 。 

另 一 方面 ， 从 频 域 的 角度 看 ， 上 述 所 列 各 类 客观 现象 频谱 呈现 一 种 と 規律 
ル / 过 程 是 一 种 非 平稳 过 程 ， 是 通过 带 通 滤波 后 表现 出 来 的 统计 自 相似 过 程 ， 因 其 
功率 谱 具 有 如 下 特点 而 得 名 。 


S (m) y -2H «1 


T 
lol’ 
.. G.Wornell 先后 在 1990 年 至 1995 SERA] Az SOS i Ae ON ay Yit 
程 的 小 波 表示 和 人 合成、 小波 系数 相关 性 等 问题 做 了 大 量 的 研究 工作 ， 得 出 小 波 变 换 
对 自 相 似 过 程 具有 好 强 的 去 相关 作用 结论 : 小 波 系数 的 相关 性 系数 随 着 采样 间 隅 增 
大 而 衰减 的 速度 要 远 快 于 分 形 函 数 本 身 目 相关 国 数 随 采 样 间 隅 增 大 而 衰减 的 速度 。 
除 此 而 外 ， 在 分 形 其 他 方面 ， 还 有 一 些 涉及 小 波 的 比较 零散 的 研究 。 壁 如 ， 小 
波 分 形 编码 四、Mallat 和 Hwang 有 关 分 形 信 号 奇异 性 的 小 波 极 大 值 方 法 等 。 
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总 的 说 来 ， 由 于 小 波 从 诞生 到 现在 不 过 十 儿 年 的 时 间 , 其 本 身 经 历 了 一 个 不 断 友 
展 和 完善 的 过 程 ， 是 现代 调和 分 析 的 一 个 重要 发 展 和 突破 ， 己 经 有 了 坚实 的 理论 基 
础 和 广 表 的 应 用 土壤 ， 已 经 为 人 们 所 接受 。 而 分 形 理论 和 小 波 相 比 虽然 起 步 比 较 
早 ， 但 因 其 “ 离 经 产道 ”的 反常 性 和 独立 性 不 易 为 大 多 数 的 人 们 所 接受 ,况且 本 身 
的 理论 基础 尚 建立 得 不 够 牢固 ， 发 展 研究 因此 就 相对 显得 碎 步 慢 行 。 小 波 分 析 工 具 
开始 进入 分 形 研 究 领 域 至 今 则 短 短 不 到 十 年 的 时 间 ， 取 得 了 一 些 可 喜 的 成 果 ， 但 分 
形 理论 中 有 许多 存在 的 问题 ， 还 有 待 借助 包括 小 波 分 析 工 具 在 内 的 手段 和 方法 去 解 
决 。 可 以 确信 ， 小 波 的 引入 犹如 催化 剂 增添 了 分 形 研 究 的 青春 活力 ， 这 是 黎明 到 来 
的 一 丝 上 曙光。 基于 这 一 理念 和 思路 ， 本 文 利用 小 波 分 析 工 具 就 分 形 理论 有 关 问 题 进 
行 了 探讨 ， 取 得 了 一 些 比较 满意 的 结果 。 


1.3 本 文 的 主要 工作 和 创新 


1. 分 析 了 分 形 曲 线 维 数 的 传统 计算 方法 及 其 缺陷 ， 提 出 了 多 分 维 分 形 曲线 的 模 
x. 基于 小 波 闹 值 去 噪 方法 ， 给 出 曲线 上 点 的 局 部 分 形 维 数 的 计算 公式 及 算法 。 仿 
真 计 算 结果 及 图 示 均 表明 该 法 是 有 效 的 。 

自从 1872 年 Weierstrass 引入 特殊 的 Weierstrass 函数 以 来 , 到 现在 人 们 有 了 分 形 
概念 后 已 经 发 现 这 类 无 处 可 微 且 连续 的 函数 是 广泛 存在 的 。 后 来 Hardy 和 Kaplan 分 
RIZE 1916 年 和 1984 年 先后 就 Weierstrass 函数 的 〈 盒 ) 维 数 问 题 展 开 了 研究 ， 探讨 
各 类 分 形 曲 线 维 数 的 计算 一 直 是 分 形 几 何 学 的 研究 热点 之 一 。 目 前 有 关 分 形 曲线 维 
数 的 计算 主要 有 量规 方法 和 周 长 面 积 法 两 种 常用 方法 ， 它们 共同 的 特点 是 ， 用 代表 
分 形 曲 线 整体 、 平 均 的 单一 维 数 来 定量 地 刻画 曲线 的 分 形 特征 。 本 人 认为 ， 其 主要 
缺点 体现 在 以 下 两 方面 : 一 是 在 忽略 细微 结构 的 不 同 的 情况 下 ， 计 算出 的 分 维 数 客 
观 上 是 一 种 比较 粗略 的 近似 结果 。Hausdorff 维 数 对 不 同 的 测量 尺码 6 是 一 个 较 敏 感 
的 参数 ， 稳定 性 的 丢失 意味 痢 维 数 参 量 原 本 价值 也 驶 无 法 体现 出 来 ， 二 是 实际 测量 
的 数据 要 受到 各 种 各 样 的 污染 ， 由 此 会 增 大 或 者 减 小 维 数 的 估计 值 ， 同 时 在 光滑 去 
品 的 人 为 操作 下 很 显然 也 会 诚 小 维 数 估 值 。 为 此 ， 我 们 有 必要 引入 多 分 维 曲 线 模式 
有 反映 分 形 曲线 的 局 部 分 维特 征 ， 这 也 就 意味 着 曲线 上 有 了 “点 的 维 数 ”概念 。 在 定 
量 刻画 分 维 方面 ， 我 们 采用 小 波 阔 值 去 噪 方法 《〈 非 光滑 方法 ) 把 分 形 函 数 从 噪音 数 
据 中 分 离 出 来 ， 通过 收敛 性 的 理论 证 明 ， 可 以 有 效 地 保留 其 分 形 维 数 ， 仿 真 算 例 结 
果 也 验证 了 方法 的 有 效 性 。 同 时 ， 我 们 也 得 到 了 基于 小 波 去 品 处 理 后 的 分 形 曲 线 维 
数 计算 式 ， 并 可 很 自然 地 应 用 于 局 部 化 维 数 的 计算 。 
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(>) E x, 处 的 U(x,,5) 邻 域内 的 局 部 分 维 数 dim( f, x ) 满 足 : 
24 1 og?” " 2) — dim(f,x,). (n>a) 

详 见 第 二 章 。 

2. 基于 小 波 变 换 在 零点 具有 的 的 保 相 似 形 特征 ， 在 已 知 分 形 插 值 函 数 图 象 情况 
F, 从 有 限 数目 个 零点 出 发 ， 给 出 寻找 一 组 仿 射 压缩 变换 使 该 函数 图 象 成 为 未 知 迭 
代 函 数 系 统 吸 引子 的 一 种 构造 方法 。 

分 形 反 问题 由 于 在 图 象 压缩 技术 中 有 着 重要 的 作用 而 成 为 目前 分 形 研究 中 的 前 
沿 课题 之 一 。 其 核心 内 容 是 寻找 由 仿 射 变换 组 成 的 迭代 函数 系统 ， 当 前 求解 反问 题 
方法 依据 包括 有 拼 贴 定理 、 和 矩 量 法 和 小 波 变 换 三 种 。 和 后 一 种 相 比 ， 前 两 种 都 涉及 
复杂 的 方程 优化 问题 ， 运算 量 大 且 稳 定性 差 。 对 于 分 形 反 问题 的 小 波 分 析 ， 可 归结 
为 一 类 特征 点 尺度 关系 的 确定 。 特 征 点 的 选择 依据 是 : ”一 是 要 能 够 完整 地 描述 整个 
分 形 内 列 的 不 变性 ; 二 是 计算 量 应 尽 可 能 小 。 Mallat Æ 1992 年 证 明了 大 多 数 函数 
的 小 波 变换 模 极 大 值 表 示 可 以 完整 表达 ， 计算 量 的 减少 是 明显 的 。 基于 这 一 结论 ， 
A.Arneodo 等 人 1994 年 选择 了 小 波 变 换 初 始 极 大 值 分 倪 点 作为 特征 点 ， Z.R.Struzik 


点 为 特征 点 。 后 者 和 前 两 人 的 方法 相 比 更 具有 计算 量 小 且 可 在 较 大 尺度 下 对 点 做 出 
判断 的 优点 。 我 们 的 方法 则 是 选择 函数 小 波 变换 零点 作为 反 演 时 的 特征 点 ， 处 理 的 
对 象 简化 为 分 形 插值 函数 图 形 ， 立 论 的 依据 有 以 下 三 条 : 

(1) Loganl"lzr 1977 年 就 曾 证 明了 ， 对 窄带 信号 而 言 ， 用 信号 本 身 的 过 零点 
可 以 完整 地 表征 信号 。 小 波 本 来 就 是 带 通 函数 。 由 定理 可 知 ， 此 时 通过 
WT,,x(!) 的 所 有 零点 途径 ,必然 可 以 恢复 到 WT,, xf)。 又 因为 小 波 变换 是 可 
逆 的 ， 这 样 必然 达到 恢复 x(1) 的 目的 。 

(2) 一 个 可 微小 波 w 的 导 函 数 y 依然 是 一 个 小 波 ， 且 有 等 式 成 立 : 


| owl ly [f a,b) 
Ob a r 
W, [HORACABREW, [/)ORESA, WHE, WEARER 


以 前 选 法 依然 能 够 获得 一 样 平等 的 地 位 。 
(3〉Barmnsleyt 叶 也 已 证 明 过 ， 分 形 插值 函数 是 定义 于 紧 区 间 上 的 自 仿 射 函 数 ， 


是 迭代 函数 系 的 唯一 吸引 子 ， 是 产生 分 形 图 形 的 重要 构造 方法 。 
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文中 证 明了 小 波 变换 零点 处 的 保 相 似 性 特征 ， 给 出 了 分 形 插值 函数 反 演 的 解 及 
算法 ， 算 例 验证 该 算法 的 有 效 性 。 
详 见 第 三 章 。 


3， 提 出 了 分 形 时 变 维 数 的 概念 ， 将 其 引入 两 类 具有 局 部 自 相似 性 的 有 偏 随机 过 
程 中 。 借 助 小 波 基因 有 的 尺度 特性 很 适合 分 析 非 平稳 性 和 局 部 自 相 似 性 这 一 特点 ， 
给 出 Hurst 指数 的 小 波 估算 公式 及 算法 ， 证 明了 小 波 估算 值 序 () 和 真实 值 (OE 
容 的 ， 并 用 仿真 结果 加 以 验证 ， 同 时 对 上 海 股票 市 场 月 指数 变化 进行 了 实证 分 析 。 

为 了 揭示 自然 界 中 一 大 类 不 规则 形体 的 标 度 不 变性 ， 一 方面 人 们 发 明了 各 种 各 
样 的 分 形 维 数 〈 由 标 度 关系 推出 的 一 些 确定 的 量 值 ) 用 来 反映 分 形 对 象 的 自 相似 构 
造 规律 ， 另 一 方面 人 们 在 研究 自 相 似 性 时 又 总 是 假定 维 数 是 一 个 常数 ， 忽 略 了 分 形 
体 在 一 维 时 间 世 界 的 演化 常会 导致 自 相似 性 的 改变 。 为 了 能 刻画 时 间 序列 分 形 的 时 
空 动力 学 演化 过 程 ， 有 必要 引入 时 变 维 数 D() 概 念 。 又 由 于 平面 空间 中 分 数 布朗 运 
动 (FBM) 中 的 Hurst 指数 和 相关 分 维 数 存在 着 D = 2-HIR-BRBRKA, ARAN we 
选取 分 数 布朗 运动 和 分 整 自 回归 滑动 平均 模型 (ARFIM) 作 为 突破 口 ， 把 传统 的 两 类 
模型 拓展 为 下 面 两 式 ， 

r0 f he-u% -uA paseo) [i-u ano) 


(i- o(8)ft - 8)" x(t) = (-e(B)el) 


此 时 D()-22-H(). 4A()= V, 时， 就 是 一 般 的 FBM 和 ARFIM 两 种 模型 。 

对 于 上 述 具 有 局 部 自 相似 性 的 随机 过 程 , 我 们 利用 Daubechies 小波 建立 了 随 机 変 
量 连续 小 波 变换 平方 与 小 波 变换 尺度 之 间 的 一 个 局 部 对 数 线性 关系 式 ， 通过 最 小 二 
乘法 得 到 HC) 的 估计 玉 (t)， 并 在 数学 上 得 到 有 关 #() 和 HH() 的 两 个 极限 表达 式 ( 相 
容 性 结果 ) : 

EBO H(). — ValB()-0. (L0) 

提出 了 实现 这 一 过 程 的 算法 ， 算 例 结果 也 显示 该 算法 是 有 效 的 。 另 外 ， 对 上 海 
股市 的 实证 分 析 ， 对 于 投资 策略 决策 有 着 重要 的 指导 意义 。 同 时 ， 我 们 相信 ， 在 目 
前 金融 工程 的 未 来 研究 发 展 中 ， 这 一 概念 及 其 小 波 分 机 方法 的 价值 将 会 愈 来 愈 快 地 
体现 出 来 。 

详 见 第 四 章 。 
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4. 借助 小 波 变换 对 自 相 似 过 程 所 起 的 白化 作用 和 最 小 二 乘法 ， 建 立 了 无 须知 道 
区 过 程 中 参数 分 布 的 谱 参 数 估计 方法 。 

区 过 程 在 生产 与 生活 中 是 一 类 很 常见 的 统计 自 相似 随机 过 程 。 其 分 维 谱 特 性 满 
ETA: 

MOLD Yy ; y=2H +1 

以 前 人 们 在 时 间 域 中 采用 ARMA 和 FBM 两 个 模型 基本 上 均 能 很 好 地 表达 这 一 
随机 过 程 ， 但 其 处 理 方法 较为 复杂 My 的 选择 也 不 近 人 意 。 由 于 该 过 程 表现 出 
的 非 平 稳 性 ，G.Worneliu5 在 1992 年 首次 给 出 了 刻画 y 过 程 的 小 波 模 型 ， 并 提出 了 
用 极 大 似 然 估 计 方 法 估算 过程 谱 参 数 y 的 途径 。 但 这 种 方法 要 求人 们 知道 该 过 程 
的 概率 分 布 函数 ， 应 用 起 来 会 带 来 有 关 分 布 的 假设 成 份 。 为 了 客观 地 估计 参数 y 
值 ， 通 过 对 多 过 程 的 小 波 变换 系数 之 间 在 相关 结构 上 的 特点 以 及 小 波 消失 矩 对 系数 
相关 性 衰减 性 的 影响 两 方面 的 研究 ， 我 们 给 出 了 基于 小 波 系数 方差 结构 的 谱 参数 / 
估算 最 小 二 乘法 ， 并 用 仿真 对 方法 的 有 效 性 进行 验证 ， 取 得 较为 满意 的 结 末 。 
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5， 提 出 在 最 优 尺 度 上 跟踪 小 波 模 极 大 值 线 的 参数 方法 ， 为 研究 多 分 形 谱 厅 异 性 
提供 了 一 种 途径 。 探 讨 了 分 形 测度 小 波 变换 极 大 值 线 的 拓扑 分 岔 情 况 ， 并 提出 了 一 
种 基于 小 波 变换 模 极 大 值 的 多 重 分 形 谱 估算 方法 。 

多 重 分 形 是 是 定义 在 分 形 上 的 、 由 多 个 标 度 指 数 的 奇异 测度 所 组 成 的 集合 ， 主 
要 研究 分 形 测 度 在 分 形 支 集 上 的 分 布 情况 一 一 多 分 形 分 维 谱 。 要 完整 地 刻画 有 关 局 
部 分 维 c 在 空间 分 布 的 信息 f(g)， 小 波 方法 是 分 析 多 分 形 局 域 结构 的 有 效 工 具 。 
Mallat 和 Hwang 9 1992 年 开创 了 利用 小 波 模 极 大 值 来 分 析 分 形 信号 局 部 奇异 性 行为 
的 先河 ， 该 法 中 的 局 部 极 大 值 点 〈 在 某 一 尺度 级 水 平 上 形成 极 大 值 线 ) BEAN 
度 上 进行 计算 的 ， 要 涉及 大 量 的 数值 运算 及 拓扑 问题 ， 特 别 是 对 于 分 形 中 非 孤 立 的 
突变 点 ， 在 跟踪 时 就 会 不 知 所 措 。 在 研究 多 分 形 谱 的 过 程 中 ， 跟 踪 失 效 也 就 意味 着 
分 形 局 部 细节 信息 的 丢失 ， 不 利于 谱 维 数 计算 结果 的 真实 性 ， 这 是 分 形 信号 与 一 般 
言 号 的 不 同 之 处 。 

本 文 假定 小 波 变 换 是 一 个 连续 变换 ， 在 一 个 最 好 尺度 上 所 有 极 大 值 点 都 可 以 得 
到 确定 ， 而 每 一 个 极 大 值 点 均 可 以 看 成 一 条 极 大 值 线 的 起 始点 。 为 此 ， 文 中 选择 上 其 
有 nn 阶 消失 矩 及 良好 局 域 性 的 Gauss 函数 作为 母 小 波 ， 把 小 波 变换 中 的 两 个 变换 因 
子 a 和 4b 看 成 是 参数 1 的 函数 alf)、2(1)， 这 样 沿 某 条 极 大 值 线 极 大 值 点 运动 的 轨迹 
可 由 下 面 的 一 个 偏 微分 方程 决定 : 


= 分 形 理 论 中 若干 问题 的 小 波 解 法 


A LT Want =0 
在 FV" =0 的 前 提 下 , 理 疹 上 可 用 Runge Kutta 方法 求解 上 述 方程 。 

在 有 极 大 值 线 分 贫 点 或 极 大 值 点 邻近 比较 平滑 的 情况 下 ， Mallat 方法 及 上 述 的 
参数 跟踪 算法 就 没有 明显 的 效果 。 文 中 选取 了 小 波 变 换 在 分 命 点 邻 域内 Taylor 级 数 
展开 中 的 两 个 多 项 式 函 数 作为 GS 范式 中 的 芽 函 数 ， 探 讨 了 这 类 奇异 点 附近 极 大 值 
ZR BE In] TE DG o 

小 波 神 极 大 值 线 结构 包含 了 多 分 形 分 形 测度 奇异 性 的 必要 信息 ， 但 由 于 分 形 的 
奇异 点 是 非 弧 立 的 ， 要 计算 每 一 点 态 奇异 指数 c(xo) 是 很 困难 的 。 文中 提出 了 一 个 
基于 小 波 变换 横 极 大 值 的 广义 g 次 信息 量 维 数 概念 ， 然 后 通过 多 重 分 形 语言 转化 式 
完成 对 多 重 分 形 谱 的 计算 。 这 一 过 程 中 ， 对 于 小 于 尺度 a 的 所 有 尺度 下 每 一 个 极 大 
值 点 处 的 小 波 变 换 模 极 大 值 采 用 沿 极 大 值 线 的 上 确 界 代替 ， 这 是 解决 非 弧 立 点 附近 
奇 寞 性 度量 问题 的 关键 步 又 。 

详 见 第 六 章 ， 
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第 二 章 多 分 维 分 形 有 曲线 维 数 计算 的 小 波 方法 
2.1 引言 


分 形 曲 线 是 自前 分 形 几 何 学 中 的 研究 热点 之 一 ， 其 研究 成 果 也 相当 丰富 。 有 
关 分 形 曲 线 的 研究 可 以 追溯 到 Weierstrass1872 年 当时 引入 的 无 处 可 微 连 续 函 数 : 


W, (x)= y cogs” nx] 0<Ai<i 
n=] 


他 曾 证 明了 当 s EPR, な >1+ ニ ァ ヶ , WIER [01], W(x) RAR 
存在 《有限 或 无 限 ) 的 导数 。 而 1916 年 Hardy WX uEB] T, 3 As21, WEE 
fs] — AW, ,(x) 均 没有 有 限 的 导数 。 在 分 形 概念 提出 以 后 ， 其 实 我 们 可 以 看 到 无 处 
可 微 连续 函数 是 广泛 存在 的 ， 于 是 人 们 转 而 将 多, ,(x) 函数 作为 一 个 具有 分 形 特征 
的 函数 加 以 研究 ， 集 中 探讨 了 函数 图 象 的 维 数 问题 。 

3S EA T(W) RR WW, ,(x) 的 图 象 ， 在 Mg >1 的 条件 下 , Hardy 在 1916 年 证 明了 
其 盒子 维 数 满足 的 一 个 结论 ; 
logA 


dim (W )x24 
ogs 


随后 1984 年 Kapian 证 明了 : 
log A 


dimT(W)-24 
log s 


设 分 形 曲线 C 是 连续 映射 SHR. fF:19"'. ARRIE SEAT 
算 ， 现 阶段 主要 有 以 下 两 种 方法 。 

第 一 种 方法 : 量规 方法 

设 f 是 连续 的 双 射 ， 显 然 C 是 一 个 无 自 交 点 的 Jordan 曲线 ， 分 形 曲线 的 各 类 
长 度 公 式 具 有 下 列 关 系 式 : 
| Los sg 
L; =Ly ig 
L ie 


由 此 可 以 推出 ， p=1-'8 P -— 2.1.1) 


HPL 表示 CC 的 欧 几 里 德 长 度 ，ZLos 表示 7 的 欧 几 里 德 距 离 ， Zr 表示 C 的 
Hausdorff 长 度 ，D 是 曲线 C 的 维 数 ，& 方 株 度 , 8 为 码 尺 ， 而 且 有 : 
=E 5 ts (2,1,2) 


— 
OC 


将 上 式 代 入 (2.1.1) 式 ， 可 得 到 下 面 的 维 数 计算 公式 ; 


des), 、 | | 
D =A ond) 213) 
AL. = Lz» VLA REPAY ALBERS S 的 两 规 脚 所 测 曲线 C 的 长 度 。 这 样 


Lif 实际 上 就 是 以 5 为 步 长 测 得 的 步 数 N(5)。 
(2.1.3) 式 就 可 以 改写 为 ， 
log N(ó 
icd y. 


jc — HC EE. HE A S PMA BE HE log 6, log") | 所 得 直线 
斜率 就 是 D。Okubo 等 人 的 采用 这 一 方法 在 6 =1-15 之 间 测 得 圣 安 德 烈 斯 断层 分 
维 D=1.31+0.02。 中 国学 者 李 后 强 等 人 叫 利 用 此 法 计算 几 种 肌 红 绰 白 蛋 白质 分 
THER Hausdorff 维 数 ， 取 得 了 一 定 成 果 。 但 这 种 方法 难免 忽略 了 蛋白 质 内 部 结 
构 的 细节 问题 ， 只 能 是 提供 定性 参考 数据 。 

第 二 种 方法 ， 周 长 面积 法 

对 于 封闭 曲线 C ， 可 以 用 C 的 周 长 和 C 所 围 区 域 面积 4 来 计算 分 维 D， 亦 岂 
小高 法 。 

*5 P. EM C 的 欧 几 里 德 长 度 时 ， 有 下 式 成 立 : 

P,” m g の の AME eme (2,1,4) 

其 中 5 RGR, a 为 无 量 刚 的 常数 。 

设 P, 为 曲线 C 的 Hausdorff KE, WA: Pp,” = a, AA 

由 (2,1,4) 式 可 以 推出 维 数 的 计算 公式 ; 


de^ o qq 0 
2 log a, + togt } (2-1,5) 
若 研究 对 象 是 一 个 几何 图 形 ， 取 不 同 码 尺 5 (i =1,2,…,n) 测 定 曲线 C 周 长 值 
A” P — 
な 和 面积 值 4 ， 然后 对 | E 73i 进行 线性 回归 即 可 求 出 刀 值 。 在 


这 里 w 和 万 均 假定 是 常数 。 此 法 在 金属 断裂 韧性 的 研究 中 发 挥 了 作用 [1。 


若 要 计算 天 个 几何 图 形 周 界 整体 的 维 数 ， 则 要 对 每 个 几何 图 形 按 上 述 方法 调 
量 ， 数 据 按照 (2,1,5) 式 回归 出 D,， 取 个 图 形 分 维 DD 的 算术 平均 值 作为 整体 维 数 
D. RA: 
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综 上 所 述 ， 对 于 一 条 分 形 曲线 C ， 用 代表 曲线 C 整体 、 平 均 的 单一 维 数 来 刻 
而 曲线 的 分 形 特征 是 上 面 两 种 方法 的 特点 。 的 确 ， 维 数 的 本 质 是 集合 层次 的 一 种 
量 值 标 号 ， 但 就 实际 中 的 蛋白 质 分 子 链 而 言 ， 由 于 其 分 子 结构 的 复杂 性 ， 要 精确 
地 刻画 其 分 形 特 征 ， 人 迫切 需要 研究 其 局 域 分 维 乙 。 在 这 种 情况 下 ， 以 往 的 计算 方 
法 就 显得 不 太 适 应 ， 其 缺点 主要 体现 在 : 一 是 在 忽略 细微 结构 的 不 同 的 情况 下 计 
算出 的 分 维 数 客观 上 是 一 种 比较 粗略 的 近似 结果 。 不 同 的 尺码 5 ， 计 算得 到 的 分 
维 数 是 一 个 较 敏 感 的 参数 ， 稳 定性 的 丢失 则 意味 维 数 参 量 原 本 价值 无 法 体现 出 
来 。 几 类 蛋白质 分 子 链 的 分 维 测定 结果 表明 了 这 一 点 中， 二 是 实际 中 测量 的 数据 
是 受到 各 种 污染 的 数据 ， 这 种 污染 “噪声 ”会 增 大 或 减 小 分 形 维 数 的 估计 值 外 |， 
并 且 对 噪音 数据 的 光滑 也 会 减 小 分 形 维 数 。 基 于 以 上 两 方面 的 原因 ， 我 们 需要 一 
个 多 分 维 分 形 曲 线 的 模型 及 一 种 新 的 计算 方法 ， 应 该 具有 解决 上 述 所 存在 问题 的 
两 个 优点 : 其 一 是 它 能 够 很 好 地 反映 分 形 曲 线 的 局 部 分 维特 征 ， 即 意味 痢 由 点 到 
点 的 “分 形 指数 ”各 不 相同 ， 其 二 它 不 仅 能 够 对 测量 数据 进行 去 品 ， 把 分 形 函 数 
从 加 噪音 数据 中 分 离 出 来 ， 而 且 能 够 有 效 地 保留 其 分 形 维 数 ， 尽 可 能 地 克服 平滑 
性 。 显 然 ， 带 有 平 请 性 质 的 去 噪 方法 是 不 能 用 的 。 

为 了 解决 这 个 多 分 维 分 形 曲 线 的 维 数 问 题 ， 我 们 引入 基于 小 波 加 阔 去 品 的 方 
法 。 

本 章 主 要 工作 安排 如 下 : 在 第 二 节 和 第 三 节 ， 对 所 涉及 的 盒子 维 数 及 相关 小 
波 方面 的 知识 加 以 回顾 。 然 后 在 第 四 节 ， 基 于 小 波 加 国 去 噪 方法， 我 们 给 出 分 形 
曲线 局 部 维 数 的 一 种 小 小 估算 方法 。 在 最 后 一 节 ， 用 仿真 例子 对 该 方法 加 以 验证 。 


2.2 盒子 维 数 的 定义 和 计算 


合 维 数 (Box Dimension) 是 应 用 最 广泛 的 维 数 之 一 ， 主 要 是 由 于 其 数学 估算 和 
经 验 估 计 相 对 其 它 分 形 维 数 来 说 容易 一 些 。 通 常情 况 下 ， 曲 线 越 粗 明 ， 邦 它 的 分 
形 维 数 就 越 大 。 

假定 G= {(x, f (x)x eoi] e 98 是 定义 于 [01] 上 的 函数 f(x) 的 图 象 。 对 于 
s>0，85(z,s) 代 表 以 z 为 中 心 、 半 径 为 的 圆 ，N_ (GO) 是 覆盖 G ND Bz.) 的 最 
小 数目 ，Z_(G) 为 E_ Minkowski 香肠 G(s) 的 Lebegue 测度 。 

G(e)= LJ a(z. £) 
zeG 


f(x) 的 分 形 维 可 以 定义 为 
dim( f) " lim log N, y," y — limf 2 _ log L{G(e ) jos ;) "M V (2, 2,1) 


在 上 式 中 N,(G) 和 ZL{G(e 上 } 直 接 分 析 较 为 困难 ， 转 而 采用 f(x) 的 変 差 来 求 其 
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分 形 维 。 
定义 fx) he- RAA: 


Osc, (x) = Suplf (,)- F(x, Jo. x; € [x —E,X+ £f [0.1] 


定义 一 变 差 为 : V(e)= | Ose, (xix 
于 是 /(z) 的 维 数 就 可 以 写成 下 式 0og， 
dim(f(x))- iml 2 ー log (e) d "LZ lul (2,2,2) 
若 (2,2,1) 和 (2,2,2) 两 式 中 的 极限 不 存在 ， 我 们 可 以 用 上 极限 、 下 极限 蔡 代 。 
另外 的 一 种 计算 途径 和 Holder 指数 有 关 。 设 是 C 一 个 正常 数 ，/(x) 满 足 指 数 
@ 的 Holder 条件 : 
If(x,)- f(x,) < Clx, 一 xil EES (2.2.3) 
结合 (2,2,2) 式 ， 可 以 推出 dim(f(x))<2-a 
另 一 方面 ， 若 存在 一 个 正 整 数 pP>1，C>0， 使 得 对 于 小 的 正 数 庆 > 0， 有 下 
式 成 立 : 


fü f(x h)- f(x - hy! ae" BO. SCs (2,2,4) 


这 也 表明 bo dim( f(x) 2 2-a. 
ER (2,2,3) 1 (2,2,4) P 5X AT DAE dim( f(x) 和 Holder 指数 联系 在 一 起 ， 并 且 确 
保 (x) 的 第 数 dim(/(x)) -2- a. 
f(x) E x, 处 的 局 域 维 数 定义 为 : 
dim(j(z)z)=limdimG 门 Boe) 
当 给 定 有 限 的 数据 集 时 ， 也 可 以 计算 维 数 。 不 妨 设 es 为 一 条 分 形 曲线 的 度量 ， 
并 定义 @(e)= "ey. CF(e) 为 曲线 的 方差 )， 让 尺度 es (E12, n) 不断 変 小 , id 


算数 据 对 集 (logs | log(Q(e, )]) 然后 求 出 回归 直线 斜率 即 为 分 形 维 数 。 
2.3 与 维 数 计算 有 关 的 小 波 分 析 


设 w(x) 是 一 个 具有 紧 支 撑 的 Daubechies 小波 
y,,(x)=2% y(2/x-k), j,keZ 


其 中 jk DANE 27 REM k 位 置 的 多 分 辨 指标 值 。 对 于 (xer (R) w 
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GE 
f (x)= 0 が (x) , C. 7 [FW (x)dx 


对 支撑 于 [0,1] 60:3 f(x), Xe FEE RE AT = 01, k 9012/71, 3f 
对 于 w,, (x) 做 一 些 限制 ， 使 得 y,, (x) 成 为 一 个 完备 区 间 正 交 基 。 

& D2S^ f(n)(ne NYA h f(x) 的 离散 采样 序列 ，W4f(n) (ne NN) 为 DD 在 
每 一 尺度 /上 的 小 波 变换 值 ，S4 7( ヵ ) (ne NN) 为 DD 在 尺度 j 上 的 近似 。 小 波 快 速算 
法 是 在 每 一 尺度 /上 把 $4 了 分 解 为 下 一 尺度 的 S24,f 和 Ws, f 

Sa SO /hs Wf = fg j501-,J-1 
其 中 了 ABIERE, ho g ERRERA olx) MR RC y(x) AM RE BO 28 


波 器 。 小 波 快速 重 构 算法 则 是 下 面 的 表达 式 : 
S4 太 = S4 fth Wi mE. je0beJ-l. een (2.32) 


基于 小 波 系数 c， 的 衰减 估计 ， 有 利于 估计 f(x) 的 光滑 度 。 
定理 2.2.1 车 一 个 有 界 函 数 f(x) 是 Holder EAN, O<asl 
I(x) FO] s clx- y 
那么 f(x) 的 小波 変換 系 数 満足 le «02799, es (2,3,1) 
证 明 如 下 : ”由 于 |w(xjzxe =0， 于 是 有 
= rtv x-ri- /lar fs 7^ v x - k]ehx -k| dx 
«C279 [yO dy 
0 
另 一 方面 ， 我 们 又 有 以 下 的 定理 (1; 着 是 y(x) 紧 支 小 波 ，/(x) ER GENUS 
数 。 对 于 某 we (01), lels cC, WA f(x) ZIM a - Holder 连续 的 。 


结合 (2,2,2) 式 ， 由 1994 年 Massopust 5! 所 推 的 结果 ， 可 知 f(x) 的 小波 系 数 可 
被 用 来 计算 曲线 的 维 数 : 
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lo 2. Ck 
dim(f(x)= ¥, + lim 4 k | リ / ees (2,3,3) 


该 结果 是 在 没有 考虑 数据 受到 “污染 ”的 情况 下 得 出 的 。 


2.4 ”基于 小 波 去 噪 的 分 形 维 数 估算 方法 


2. 4. 1 分 形 曲 线 的 小 波 加 羡 去 品 处 理 


利用 小 波 变换 进行 去 噪 及 重 构 一 直 是 一 个 热门 话题 。1992 年 Mallat WA TA 
异 信号 及 噪声 在 小 波 变换 各 尺度 空间 中 模 极 大 值 的 不 同 传播 特性 。 该 方法 经 过 对 
模 极 大 值 处 理 之 后 ， 存 在 一 个 由 剩余 模 极 大 值 重 构 小 波 系数 的 问题 ， 计 算 比 较 复 
杂 ， 对 原 有 曲线 有 一 定 平滑 作用 [9。 另 一 种 小 波 去 噪 方法 hl 是 直接 对 信号 的 
小 波 系数 采用 软 或 硬 ) 阐 值 方法 进行 “过 滤 ”， 计 算 量 小 ， 方 法 较 简 洁 ， 同 时 
也 不 会 过 分 地 改变 曲线 的 光滑 度 ， 所 以 适宜 于 处 理 分 形 曲 线 。 

假设 数据 y, yos y, 是 从 下 面 的 非 参数 回归 模型 获得 的 


= ft i212,-.2" a (2,4,1) 

其 中 x, = eu]. e, ~ NO ruo 的 高 斯 随机 噪音 ， 7(x) 是 定 叉 
于 [ojJj 上 的 未 知 待定 的 函数 。 

在 去 噪 过 程 中 获得 估 值 f(x) 的 算法 如 下 : 

步骤 1: 利用 上 节 所 述 的 Daubechies 紧 支 小 波 滤 波 算法 来 处 理 观 测 数 据 y 获 
(Bee NE y,,, j=0,1,…,N-1， 上 =0,1,…,2/ 一 1 

步骤 2: 通过 下 面 的 软 阐 规 则 对 经 验 小 波 系 数 y， 进行 加 限 处 理 得 到 
56). 

5, Je bul) al 2 t, 


m くら 
EU LE, = 9/ [- A2 logín) e 
n 
歩 3， 由 加 限 经 验 小 波 系数 5. (y,, ) 通 过 (2,3,2) 式 重 构 获得 f(x, ) 的 小 波 估 
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计 值 f(x,)。 

可 以 看 出 ， 选 取 一 个 立 值 1, 后， 把 低 于 4, 的 小 波 系数 y,。 (这 部 分 主要 由 品 
He 引起 ) 设 为 零 ， 而 高 于 4 的 小 波 系数 y,，。( 主 要 由 (x,) 引 起 的 ) 经 过 处 理 后 
的 8 (y,, )， 就 可 基本 上 理解 为 是 由 (x,) 引 起 的 ， 进 而 得 到 一 个 估计 了 (x,)。 


2. 4. 2 多 维 分 形 曲 线 局 部 维 数 的 计算 


在 (2,3,3) 式 成 立 的 条 件 下 ， 我 们 考察 f(x) 在 x 处 的 邻 域 U(xo,5) 内 的 情况 得 
到 f(x) 去 噪 后 的 一 个 局 部 化 维 数 结论 : 
4n> oft, 在 概 率 意 又 上 
p Liog 5 $5. (v, " J 5dim(f(x)x,)) ^ …… (2,4,2) 
(n) 
"Asc SP diam a a adim(f x) n] p 


某 一 正 数 ， 7 = 01 ,J,k 20,527 -1. 
证 明 如 下 : 


: | | 
Ü Vip 7C, tS j201,J, kz04,2/ -1. ……(2.4.3) 
其 中 アル sC k XP 别 是 y, f(x ) £, 的 经 验 小 波 系 数 。 

此 时 ， 加 限 小 波 系数 5, (v, Es 


> AE 3 (v,.]2 2 Yal h cen (2,4,4) 


EFIE, HEH o 的 独立 正 态 随机 变量 ， 而 小 波 变换 是 正 交 的 ， 所 以 
,也 是 零 均值 、 方 差 为 0? 的 独立 正 态 分 布 。 


Chebyshev 不等式 表明 : Se, | 为 一 个 2 次 多 项 式 
当 n ->w% 时 ， 现 在 考察 (2,4,4) 式 两 边 : 


oe Dy, il|- oe: 5k | 


= log, DM +0, (n 27^ ) 
k 


1 
C yk UAR 
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E: ee (Y... | + 08 + の py lea | 
k 
= ke (Y. + 0, 2^ [Y a | y ten (2,4,5) 
k 
同时 有 


og: Ena -2” 3 . e [e + O,(n #2h )- 2^ n} 
k k 
= log, pA | * log, | 十 o, [tos ny? n ^ 2 Fkn] 
= ee Y. | +0, C ny? n^ 2^ Xle, k | TELLE (2,4,6) 
上 面 各 式 中 的 O。 均 表示 Taylor 多 项 式 展 开 。 
log; 3765] 
再 由 (2,3,3) 式 得 到 ; 2dim(f)=3+2 lim — ! f 
ops[ le, sl | -Eam)-3+ ote, 


ye E [ism «097 l B(n) ES O(n EN oo) 
k 


于 是 有 
(log; ED C,. ; "E Uo 
pp iei) on ua = im oes 1°) ETIN, 


Á 
= in (r 2 | っ E-29m(7 た (2 | 
n» Jn 
P Ny 
APD Ju 


将 上 式 代入 (2,4.53)、(2,4,6) 两 式 ， 并 利用 (2,3,2) 式 结果 我 们 可 以 得 到 下 式 ; 
og: っ dm の )- テ . (n2) 
log. Psal 2^1. i  dim)-2. (n>) 

由 (2,4,4) 不 等 式 关 系 ， 我 们 就 可 得 到 和 (2,4,2) 式 类 似 的 维 数 计算 公式 ; 


A SIR eM NOT B 


oo 


: * log; (v, 5, (v, A x dim(f, Xa) j (n =y oo) 


以 上 结果 可 以 说 明 以 下 几 点 ; 
(D 在 无 噪音 情况 下 ， 在 高 分 辨 级 水 平 上 ， 小 波 系数 可 以 刻画 f(x) 的 局部 
细节 特征 ， 从 而 确定 函数 /(x) 的 维 数 。 一 旦 随机 噪音 进入 ， 则 它 在 更 
高 分 辨 级 上 左右 小 波 系数 c，,， 使 得 ck 失真。 对 于 满足 指数 w - Holder 
条件 的 f(x), ， 此 时 ， 当 (2.3,4) 式 中 的 噪音 项 普 4e，, 达 zz“ 阶 数量 级 时 ， 
我 们 可 以 看 到 在 更 高 分 辨 级 上 ，y,, 主要 由 a,, 控制 ， 分 形 曲线 f(x) 更 
详细 的 特征 也 就 难以 很 好 地 恢复 起 来 。 

在 相对 较 低 的 分 辨 级 j KEE, y WEZH ,控制 。 正 象 上 述 定 
理 所 证 ， 对 于 所 取 较 大 的 样本 数 n，/(x) 的 分 形 维 数 都 几乎 保留 在 分 
辨 级 六 的 小 波 系数 5 (y, ) 中 。 

通常 情况 下 ， 采 用 光滑 方法 处 理 曲线 (x) 后 的 维 数 估计 值 为 1， 也 就 
是 说 ， 用 光滑 办 法 恢复 原 函 数 的 维 数 是 不 可 能 的 。 

这 种 计算 方法 既 可 用 于 计算 f(x) 的 整个 曲线 维 数 ， 也 可 以 在 x 的 邻 域 
U(x。,6) 内 进行 样本 点 的 采集 ， 按 照 以 上 的 算法 及 定理 ， 从 而 得 到 (x) 
在 x 处 的 局 部 维 数 估计 值 。 另 外 ， 对 于 加 阐 处 理 后 的 复原 f(x, ) 也 可 直 
接 采 用 第 二 节 末 所 述 方法 进行 维 数 计算 。 


Ne 


(2 


Ad 


(3 


Næ” 


(4 


25 数值 仿真 算 例 


算 例 1 选取 Gauss 随机 过 程 中 的 一 个 样本 轨道 曲线 作为 分 形 曲 线 f(x). B 
设 随机 过 程 用 下 式 表达 : 
X(x)=X, + fe = の の dB,, xelo] 


其 中 X, 为 一 个 标准 正 态 随机 变量 ，8B, 为 一 个 标准 的 Brown 运动 。 

” 随 着 曲线 上 点 变化 而 变化 的 自 相 似 参 量 H(x)〈 其 作用 等 间 于 维 数 ) 设 定 为 一 
线性 函数 : 厅 (x)=0.2+0.7x， 样 本 数 n = 2100 。 

图 2.1 所 示 为 /(x) 的 多 分 维 分 形 曲线 ， 图 2.2 则 为 (x) 的 局 部 维 数 变换 。 


w -% k 
算 例 2 ”对 于 Weierstrass MM: f(x)= $2) al (2) J 
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该 函数 具有 自 仿 射 性 和 相同 的 局 部 维 数 dim(D=1.5， 所 添加 噪音 为 高 斯 白 品 
声 ， 均 匀 取 值 的 样本 数 n = 800 ， 选 具有 2 Erit AN Daubechies 小 波 作为 处 理工 
具 。 结 果 见 以 下 图 2.3~2.6。 

图 2.3 和 图 2.5 结果 表明 ， 加 闭 去 品 法 对 原 曲 线 的 不 规则 程度 影响 不 大 ， 所 
测 维 数值 (1.47) 和 真实 值 (1.50) 很 接近 ， 结 果 比 较 理想 。 而 图 2.3 和 图 2.6 表明 , 
采用 小 波 模 极 大 值 方法 去 品 ， 盒 维 数 结果 (25) 则 和 真实 值 偏差 较 大 ， 说 明 有 平 
滑 作 用 的 滤波 去 噪 方法 对 原 分 形 曲线 的 恢复 不 太 合 适 。 图 2.4 显示 加 有 噪音 的 数 
据 样 本 计算 出 的 维 数 会 明显 地 改变 原来 曲线 的 维 数 。 


2.66 本章 小 结 


在 这 一 章 中 ， 我 们 首先 分 析 了 分 形 曲线 维 数 的 传统 计算 方法 及 其 两 个 缺陷 : 
一 是 忽略 曲线 局 部 细节 用 整体 、 平 均 的 单一 维 数 表达 曲线 的 不 规则 程度 ， 二 是 实 
际 测 量 数据 受到 “噪音 ”污染 后 会 增 大 或 减少 维 数 佑 值 。 针 对 以 上 两 个 问题 ， 文 
中 提出 了 一 个 多 分 维 分 形 曲线 模式 ， 强 调 了 曲线 局 部 维 数 的 存在 ， 以 便 能 完整 准 
确 地 刻画 分 形 曲 线 的 不 规则 性 。 其 次 ， 为 了 把 分 形 函 数 从 噪音 中 分 离 出 来 并 尽 大 
可 能 地 保留 其 维 数 ， 文 中 采用 了 加 国 小 六 去 品 方 法 对 样本 数据 进行 处 理 ， 给 出 了 
利用 小 波 系 数 计 算 曲 线 的 局 部 分 维 数 公 式 ， 并 证 明了 其 收敛 性 。 

在 文章 的 最 后 ， 用 两 个 仿真 算 例 对 于 以 上 所 提问 题 加 以 说 明和 验证 。 


—— 
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Y(t) 
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9.0 | 
E21 f(x) 的 多 分 维 分 形 曲 线 ， 


H(t) 


0.6 0.8 1.0 


0.2 0.4 
图 2.2 f(x) 的 局 部 维 数 变换 。 


0.0 
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K 2.3 原 weierstrass 分 形 曲 线 (djm(D=1.5) 


fi» 


图 2.4 MARA JMA BE dim(f)-1.96) 


f(t) 


fi» 
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图 2.5 IRA ee Cdim(f=1.47) 


图 2.6 模 极 大 值 去 噪 后 的 曲线 (dim( 们 =1.24) 
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三 章 ”分形 反 演 问题 的 小 波 解法 


分 形 反 演 问 题 是 目前 分 形 研究 中 的 前 沿 课题 之 一 ， 也 是 最 为 困难 的 问题 。 分 
形 反 问题 就 是 找 分 形体 所 具有 的 各 种 不 变性 ， 这 些 不 变性 常常 和 一 些 数学 变换 庄 
如 缩放 、 平 移 和 旋转 等 相 联系 ， 可 以 说 解决 反问 题 就 是 等 价 于 找 出 这 些 基 本 的 数 
学 变换 来 。 而 一 个 分 形 集 的 产生 往往 是 在 上 述 不 变 式 基 础 上 通过 友 代 过 程 来 实 
现 ， 所 以 迭代 函数 系统 是 构成 分 形 反 问题 的 重要 部 分 。 


3. 1. 1 分 形 反 问题 研究 发 展 的 历史 概要 及 主要 应 用 


在 1985 年 M.F,Barnsley 中 利用 一 类 设 定 的 迭代 函数 系统 来 通 近 分 形 的 不 变 
性 ， 从 而 开创 了 分 形 反 演 问 题 研究 的 先河 ， 但 方法 并 不 是 对 分 形 本 喘 所 具有 真实 
不 变性 进行 分 析 只 能 认为 是 分 形 反 演 工作 的 一 次 尝试 。A.Arneodo 等 人 在 1994 年 
至 1995 年 期 间 利用 小 波 变换 和 拼 贴 定理 在 对 一 维 “ 饼 格 ”(Cookie-Cutters) 和 和 角 问 
集 DLA 的 反 演 中 获得 成 功 趾 ， 实 现 了 这 一 领域 一 次 真正 的 突破 ， 从 而 葛 定 了 分 
形 反 演 的 理论 基础 。 在 实际 应 用 中 ， 他 们 把 小 波 变 换 最 初 极 大 值 分 倪 点 作为 特征 
点 从 使 反 演 过 程 具 有 了 可 操作 性 。 在 这 里 A.Arneodo 等 人 并 未 对 整个 分 形 反 问题 
作 全 面 的 分 析 和 考察 , 另外 这 种 特征 点 的 选择 也 是 有 探讨 的 余地 可 言 。 随 后 Forte.B 
等 人 卜 几 乎 同期 则 用 迭代 函数 系统 来 解决 逼近 函数 图 象 的 反问 题 ， 这 种 逼近 方法 
EH) BRR RET ER. HPI RH RA 
同 ， 即 使 是 一 个 简单 吸引 子 ， 上 述 方法 也 不 会 找到 实际 正确 数目 的 映射 。1995 年 
Struzik.Z.R. 的 提出 了 利用 小 波 变换 分 析 自 仿 射 分 形 反 演 问题 的 理论 基础 ， 但 他 只 
是 给 出 示范 性 表述 ， 并 未 给 出 严格 证 明 。 另 外 ， 值 得 一 提 的 是 ，1997 年 中 国学 者 
朱 治 军 和 李 后 强 俩 人 和 3 则 从 小 波 变 换 最 末 极 大 值 入 手 ， 对 解决 分 形 反 演 问题 进行 
了 一 次 成 功 的 尝试 。 

分 形 反问 题 研究 在 图 像 和 数据 压缩 技术 方面 有 着 非常 重要 的 应 用 ， 越 来 越 受 
到 各 界 人 士 的 高 度 重 视 。 假 如 可 以 对 给 定 吸引 子 加 以 计算 ， 则 用 有 限 个 少量 的 压 
缩 变换 可 以 决定 高 度 复杂 的 精细 结构 ， 这 样 复杂 的 图 象 便 可 以 有 效 地 进行 编码 储 
存 和 传输 。 整 个 图 形 集 就 可 以 用 有 限 个 映射 进行 编码 。 这 方面 ， 数 学 家 Barnsley 
和 Sloan 1988 年 在 其 文章 中 宣称 ， 利 用 他 们 的 方法 对 静态 图 像 压缩 可 获得 高 达 的 
压缩 比 ， 这 一 报道 立即 在 工业 界 引 起 极 大 的 震动 。1992 年 软件 开发 商 Microsoft 
公司 运用 分 形 技术 推出 一 套 多 媒体 百科 全 书 在 一 张 光盘 上 可 以 储存 大 量 的 文本 、 
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图 片 、 声 音 和 动画 等 数据 信息 。1996 年 3 月 Barnsley 领导 的 迭代 系统 公司 在 所 在 
地 美国 亚特兰大 市 召开 了 一 个 “分 形 图 像 压缩 技术 专题 研讨 会 ”"， 美 国 一 总 统 亲 
临 会 场 祝贺 ， 以 示 政 府 对 此 项 技术 的 高 度 关 注 。 


3.1.2 目前 分 形 反 问题 研究 的 方法 和 特点 


在 分 形 反 问题 中 大 多 考 虚 下 列 特殊 形式 的 IFS ， 
F={X:W= (w, w, Wy}, 其 中 压缩 映射 个 数 NN 固定 ,每 一 个 w :了 一 了 为 下 
面 仿 射 变换 : 


se OF seem 


令 向 量 了 = (5 ,i=12,.3,4)， 向 量 了 = (=12)， 于 是 这 N 个 仿 射 变换 {5,7) 就 
构成 1FS 全 体 。 若 考虑 概率 分 布 声 ， 则 参数 系统 (8,P, p) HORS IFS. 

目前 IFS 分 形 反 问题 的 求解 概括 起 来 有 下 面 三 种 方法 : 

(1) 基于 拼 贴 定理 的 IFS 反问 题 : 

= A 为 现实 图 像 ， 据 拼 贴 定理 ， 如 果 存 在 IFS{X;w} ， 使 得 

d, (4. | w (3) 充分 小 , R 旭 4 和 7 が 吸引 子 4, 的 距离 dw(4,4,)<(1-s) d, 也 
充分 小 。 所 以 人 们 仅仅 考虑 IFS 作用 在 现实 图 像 上 是 否 会 有 很 小 的 变换 。 

确定 型 IPl: VA, RETER, HA 

7= ala U” w, C). 
概率 型 IP2: WA, RET pjeR™, HA 
B= ad Xp (Ao w’) 

(2) Æ FEDH] IFS 反问 题 ;: 

Wt 4 为 现实 图 像 ,{G,, = | x" y" d4m» 0,n2 0] HRANIE, (g, 1» 0, » 0] 
表示 不変 測度 BI IB. 

现在 欲求 概率 IFS (X; w;p). 使 得 下 式 送 到 最小 : 

f. p)- TD" (g, (9, 3)-G,) 一 最 小 

以 上 两 类 方法 都 把 图 形 压 缩 问 题 转 化 为 求解 复杂 的 非 线性 优化 问题 ， 其 核心 
是 寻找 合适 的 优化 算法 。 璧 如 ，Barmsley 及 其 Iterated System 公司 许多 专利 产品 
都 广泛 采用 了 模拟 退火 优化 算法 。 其 大 体 上 具有 计算 量 大 、 稳 定性 差 的 特点 。 
(3) 基于 小 波 变 换 的 IFS 反问 题 

因为 实际 测量 中 涉及 的 量 是 测度 ， 故 定义 测度 u(x) 的 小 波 变 换 为 


2 分 形 理论 中 若干 问题 的 小 小 解法 


W, ula,b) = [v6 - b yu). a»0,be9i 
则 有 W, u(a,b)- Los W, u(a b) 


a =u,a ん =12,…, N. (3.1.1) 
b =u,b+y, 

WATE RS BET NEA T. AFOSR ARS, $e E CR DID 这 
样 的 对 应 点 是 相当 复杂 的 工作 。 一 般 的 解决 思路 是 : Q.1.1) OS ELI PES A 
成 立 ， 则 应 该 对 某 类 特征 点 也 成 立 。 这 样 反 问题 就 归结 为 这 类 特征 点 的 尺度 变换 
关系 的 确定 ， 特 征 点 的 选择 是 反 演 工作 中 的 一 个 重要 方面 。 选 择 依据 主要 有 两 条 : 
一 是 能 否 完整 地 找 述 整个 分 形 的 不 变性 ， 二 是 计算 量 的 大 小 。Mallat 在 1992 年 
就 已 证 明 ! 引 ， 对 于 大 多 数 函 数 来 说 ， 其 小 波 变换 的 模 极 大 值 表 示 是 f(x) 的 一 个 
完整 表示 ， 这 样 借助 极 大 值 表 示 就 能 用 较 少 的 计算 量 描 述 分 形 的 尺度 不 变性 。 有 
关 特 征 点 的 选择 ，A.Arneodo 等 人 选择 初始 极 大 值 分 岔 点 作为 特征 点 ，Z.R.Struzik 
1995 年 则 取 小 波 变 换 极 值 分 命 点 为 特征 点 ， 朱 治 军 、 李 后 强 则 取 最 末 极 大 值 点 为 
特征 点 。 后 一 种 方法 和 前 两 种 方法 相 比 具有 计算 量 小 且 可 在 较 大 尺度 下 对 点 做 出 
判断 的 优点 。 


3. 1. 3 ”本章 所 用 方法 的 思想 依据 及 工作 


本 文 所 做 的 工作 是 基于 以 下 三 个 方面 的 依据 进行 的 : 

1. Logan 曾 在 1977 年 to 证 明 ， 对 窄带 信号 而 言 ,用 信号 本 身 的 过 零点 可 以 
完整 地 表征 信号 ， 只 是 表征 不 够 稳健 而 已 。 小 波 本 来 就 是 带 通 函数 。 若 小 波 是 理 
TN, EATS, EE EWT OM MALE, LRAT LSS WT, x(t). 
X25 Bee we RE ATHY, RURARKA x(r) BO B f. US UE UE IS EAR 
不 完全 满足 ， 所 以 重建 也 就 有 点 不 够 稳健 是 可 以 想到 的 。 

2. 车 一 个 可 微小 波 y 求 导 后 w 依然 是 一 个 小 波 ， 因 为 有 下 面 的 等 式 成 立 : 


の = うみ [/ Kob) 


ap An w, [7] EO BEBCK f SREY. Lr] Re A, 也 就 是 说 ， 选 零点 为 特 
征 点 和 以 前 选 法 依然 能 够 获得 一 样 平等 的 地 位 。 

3. Barnsley 也 已 证 明 过 ， 分 形 插值 函数 是 定义 于 紧 区 间 上 的 自 仿 射 函 数 ， 是 
迭代 函数 系 的 唯一 吸引 子 〈 连 续 函 数 )， 是 产生 分 形 图形 的 重要 构造 方法 。 
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本 文 基于 小 波 变换 在 零点 的 保 相 似 性 特征 ， 在 已 知 分 形 插值 函数 图 形 的 情况 
下 ， 从 有 限 个 小 波 变换 零点 出 发 ， 给 出 一 种 寻找 一 组 仿 射 压缩 变换 
w,(i =1,2,…, 入) 使 得 该 函数 图 形成 为 待 求 JFS(W) 吸 引子 的 构造 方法 。 

在 文章 第 二 节 里 ， 考 虑 在 完备 的 Hausdorff 测度 空间 里 分 形 插值 函数 的 构 运 
问题 。 在 第 三 节 主 要 讨论 了 小 波 变换 零点 处 的 保 相 似 性 特征 。 在 随后 第 四 和 第 五 
两 节 中 ， 文 中 给 出 了 分 形 插值 函数 反 演 的 解 及 算法 并 用 算 例 对 算法 的 有 效 性 加 以 
验证 。 


32 ”分 形 插 值 函数 的 构造 


设 昌 (XY) 为 完备 空间 X 的 一 个 非 空 紧 支 集 集合 ， 给 XY 附加 一 个 Hausdoff 测 
度 , 空 同 ( な (Y)g。 ) 即 为 完备 的 测度 空间 . 
Kw X > X,i=1,2,--,N, MERIT X ERE, 即 存在 0<s <l, 
使 得 d(w,(x),w,(y))ss,d(x.y), WF x, y e X MEDR. ER MBA 


射 定理 可 知 ， 存 在 唯一 的 非 空 子 集 Bc H(X)，B = Che (B) : 


这 时 了 就 称 作 迭代 函数 系 EFSfw i 12, , N 的 吸引 子 。 
De H(R), mw =1,2,…,N): fo1]x D — [01]x D. 


"C0 LS 


Er, 2,8,y,56,0, €8,0«a, «10 «|y, «1 
4 P - (0,y,) 99 m 的 不 动 点 ,Pw m yu) Aw, 的 不 动 点 ， 选 择 参数 使 得 
w(Py)= wal) SH) i =1,2,---N 一 1 均 成 立 , 且 w(PR)w(BR)…ww(B) 不 全 共 
线 。 这 样 可 保证 [w(P)w(P,,,)| 形 成 一 非 平 凡 多 边 形 曲线 。 
另外 a ,d, 的 选择 要 确保 映射 w, 对 于 某 一 度量 来 说 是 压缩 映射 。 
Ji REX S RAR {w,,w,,…ww} 的 吸引 子 就 是 一 个 无 处 可 微 的 连续 函数 了 的 图 
形 . 函 数 f 就 是 过 插值 点 w(P)w,(P》…wv(P)》 Py 的 分 形 插值 函数 。 
对 于 一 个 函数 e [0,1], 定义 算 子 H:Cl0,1j 一 clo: 
(ur Xx) = BI (x) + y. fü Q9) 0. L6) = we») 
w(R)- P, vx elx,,,x,| 
设 是 具有 插值 点 PA A B] RR dS Ek TE wR, W 
f.) Gf. Xx) 9 f. (n 9 0) f ÈH R7 zh rx, B2) JE ET ER DICT DUBIE X — ER 
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数列 的 一 致 收敛 极限 函数 。 
另外 ， 由 于 P,P,…, Py 不 共 线 ， 则 由 JFS(w,,i=1,2,…,N) 形 成 图 像 的 分 形 维 
数 是 下 列 方程 的 唯一 解 门 ; 
3 ride?" =1 


3.3. 小 波 变换 零点 的 自 相 似 特征 讨论 


设 小 波 函 数 y 是 一 个 n+1 次 可 微 且 紧 支 于 |[-K.K| 光 滑 可 微 函 数 (n > 2) g(x) 
的 k 次 导数 ， 了 是 一 分 形 插值 函数 。 定 义 了 在 w 下 的 小 波 变 换 为 : 


Hity ,车 y 具 kk BRK 22), W= fw 为 一 迭代 函数 系统 
IFS,f£XX (b,a)e RxR,, 使 得 (-aK +b,aK -b)c [oj 成 立 ， 那 么 就 可 以 得 到 小 波 
变换 零点 的 自 相似 性 关系 方程 : 


W, f(b,a)= ーー 76 *&,ma) izl2,-.N. 


ioe tz "T Ec 


Wr ante 
2. 


P, OK +b d 2 
fairs [4 一 


W, f(b,a) = x. i /= 


] l "T +b) Fa (2 


= V 
a; Ja 1 (-aK +b) a 


l {ak +b) 


] 
Gy, Va l(-aK +b) 


A LA T (ye 


Yd P a 


鉴于 y B-TEDAI ONAN DRE, BOBS TMA: 
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we の ニー テー vf ehh e 


Oy, Ti i, (-aK +b) 


= ey. fla,b + ,,a,a) 
ay 


tf l 


于 是 得 到 如 下 命题 : 


W 
W, f(b.a)- o a) ーーー (6,q)=0 等 价 于 下 述 两 个 等 式 


W, f (a,b + €,,0,a) = 0 
SWAT 4 54 6,,0,0)= "T. | f\ab+6,,@,a)= 0 
Ob 2 r 
在 这 里 ， 定 义 A 为 满足 上 述 两 式 的 零 元 素 (ba) 的 集合 。 


w, [ke,a)=08 ee e die | 


Ay, f)= le a)e 8x9 


3.4 分形 插 值 函数 反问 题 的 解 及 算法 


WW = wow, wu) E EB SR Lo, B, y, 6, 0, 的 一 迭代 函数 系统 ， 


fA 


相对 应 的 分 形 插值 函数 ， 同 时 取 {b,a)e 4， 则 存在 点 也 ,ajes 4i=12…N. 使 得 


下 列 等 式 成立 : 
の , = Qa..……..(3.4.1) 
b=0,b+6,.…(3.4.2) 
N 
Sia, 二 (3.4.3) 
a, +E, =E (3.4.4) 
间 时 有 北 推 式 : 


(ba /)e A; AA b; = a,b" *t&,; = aa, j = 2,3,:-- 


利用 以 上 这 些 特 性 ， 可 给 一 个 分 形 插值 函数 图 形 编码 ， 从 而 得 到 最 小 数目 的 


仿 射 映射 wi Wasis Wy 
具体 算法 如 下 : 
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步骤 1: pi (b°,a°)e A, Ana? = max{a e 91, [35,(5,a) & 4} 

步骤 2; 继续 从 尺度 a 到 下 一 尺度 a? PAA (OL a! )e 4 

步 又 3， 由 计算 G.1)、(3.2) 参 数 a,s, 并 检查 是 否 符合 (3.3)、G3.4) 条 件 。 


然后 由 f 在 XX 轴 不 共 线 的 三 个 值 解 方程 ， 可 求 出 8,y,,9 的 值 


PRA: 必要 时 进行 下 一 尺度 计算 ， 重 复 以 上 步骤 。 
该 方法 的 特点 是 在 所 选 只 度 上 利用 较 少 量 的 小 波 变换 的 零点 来 构 抬 出 
分 形 插值 函数 的 迭代 函数 系 。 


3.5 ”举例 
选取 一 个 二 次 光滑 样 条 函数 作为 尺度 函数 g(x) 来 构造 小 波 函 数 g(x) 


x? O0<x<l 

-2x! - 6x -3 1<x<2 

(x)= (3 - xy 2<x<3 
0 其 它 


给 出 一 族 仿 射 映射 F = {wi,w,} 
0-5 MID 
| E " s ilt ul j im 


经 由 小 波 变换 零点 计算 出 的 参数 如 下 表 所 未 : 


|: |a] ぁ gl n. 
0.65 


Fy DUE HE SE HE HB Ra E BE E E HU TA 2 TE TS E PR CET E S LER 
数 % IFS. 
上 述 所 用 分 形 插值 消 数 、 小 波 变换 零点 及 反问 题 构 造 过 程 见 下 面 图 3.1 和 图 3.2 
所 示 。 
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3.6 AE 


在 用 小 波 处 理 分 形 反 问题 理论 方面 ， 不 同 于 以 往 特征 点 的 选择 ， 我 们 提出 了 
借助 小 波 变换 零点 这 类 特征 点 构造 分 形 插值 函数 吸引 子 IFS 参数 的 新 方法 。 首 先 
我 们 讨论 了 在 Hausdorff 测度 空间 中 分 形 插值 函 数 的 一 般 构 造 理论 ， 然 后 证 明了 
小 波 变换 零点 的 保 相 似 性 ， 给 出 一 个 零点 集 4(,) 和 一 狙 参 数 遂 推 夫 系 式 。 算 
例 表明 该 方法 是 有 效 可 行 的 ， 而 且 可 以 确定 IFS 中 映射 的 正确 数目 。 


分 形 反 演 中 存在 的 问题 及 进一步 发 展 的 方 问 : 
1. 由 小 波 变换 值 的 关系 反 推 分 形 的 变换 关系 , 然而 , 这 种 反 推 是 否 确定 叭 一 的 
分 形 ? 这 是 一 个 尚 待 解决 的 问题 。 
2， 对 于 非 相 似 分 形 ( 如 圆 和 椭圆 、 球 和 椭 球 〉 的 小 波 分 析 反 演 ， 这 必然 涉及 
尺度 参量 应 该 为 n 维 对 角 和 矩阵 的 小 波 函 数 。 
g(x)= Y ao(4Ax-b), AA MAE. 
bez’ 
这 种 小 波 函 数 的 研究 已 有 不 少 路 引 。 尽 快 地 将 其 用 于 这 类 分 形 ， 这 也 是 以 
后 反 演 工作 的 一 个 研究 方 同 。 
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40 
e 
| 
J^ 
X 
图 3.1 mw, = 1,2) BRE mt TESI Ef (x) 图 形 
200 
(b(O) a(0)) 
ば lw (bD, a (2 
(b (1), : 1) 
y * 
(b,i2,a,)) | | | | | | 
EA AAAA MARIA add MANN 


0 1 


图 32 jx) 的 小 波 变换 零点 及 由 他" ,aa yr (bY), a eite 
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第 四 章 ”时 变 分 形 维 数 的 概念 与 计算 
41 引言 


分 形 概念 揭示 了 自然 界 一 大 类 无 规则 形体 的 标 度 不 变性 ， 即 就 是 在 分 形 几何 
体 动力 学 演化 过 程 中 某 一 时 刻 在 一 定 的 标 度 尺度 范 围 内 ， 其 相应 的 测度 不 随 尺度 
的 改变 而 发 生变 化 ， 其 局 部 与 整体 在 形状 、 结 构 等 方面 具有 和 外 相似 性 。 人 们 发 明 
了 各 种 分 形 维 数 ， 它 们 是 在 分 形 意 义 上 由 标 度 关系 推出 的 一 个 确定 的 数量 值 ， 除 
了 反映 该 结构 的 自 相似 构造 规律 外 ， 尚 不 能 完全 揭示 事物 演化 的 时 空 动力 过 程 。 
与 此 同时 ， 人 们 在 研究 许多 自然 现象 的 相似 性 时 ， 又 总 是 假定 维 数 是 一 个 常数 。 
事实 上 ， 自 然 现 象 在 一 维 时 间 世 界 的 演化 常常 会 导致 自 相 似 性 变化 。 在 这 样 具 有 
局 部 自 相 似 过 程 的 情形 下 ， 上 述 维 数 不 变 的 思维 就 是 不 现实 的 。 为 了 更 充分 、 完 
整地 描述 这 种 随时 间 演 化 的 现象 ， 文 中 提出 了 新 的 概念 一 一 时 变 维 数 函数 D(x)。 

这 一 思想 受到 下 面 两 个 事实 的 启发 :彼得 斯 1991 EMA R/S 分 析 法 ( 重 标 
RER) 计算 了 分 数 布 明 运动 (FBM) 中 医 斯 特 (Hurst) 指 数 〈 反 映 序 列 参 差 不 齐 程 
度 及 其 相关 性 ) 在 不 同时 间 增 量 段 (日 ， 月 ， 年 ) 上 的 结果 ， 发 现 Hurst 指数 是 
不 同 的 这 一 事实 。 在 基因 科学 中 ，DNA 链 编码 域 中 的 核 甘酸 序列 和 非 编码 域 的 序 
列 在 自 相似 性 特征 方面 维 数 区 是 有 区 别 的 ， 若 采用 时 变 维 数 就 可 以 很 好 地 描述 
DNA 序列 在 时 空中 的 合成 演变 过 程 。 

在 平面 空间 内 ， 由 于 分 数 布朗 运动 模型 (FBM) 及 分 整 自 回归 滑动 平均 模型 
(ARFIM) 中 的 Hurst 指数 H 和 分 形 维 数 D 有 着 下 述 的 数量 关系 : D=2-H, ii 
县 分 数 布朗 运动 在 描述 地 貌 生成 、 粒 子 随机 运动 、 资 本 市 场 股 价 波动 等 方面 有 广 
泛 的 应 用 。 因 此 , 本 文选 取 FBM 模型 ARFIMA 模型 作为 突破 口 转 而 研究 时 变 Hurst 
指数 这 一 问题 。 此 时 有 D()= 2 H(t)。 同 时 ， 我 们 注意 到 ， 由 于 具有 局 部 自 相似 
随机 过 程 中 的 Hurst 指数 是 时 变 的 ，FBM 模型 及 ARFIM 模型 随机 过 程 及 其 相关 
的 增 量 过 程 一 般 来 说 都 是 非 平稳 的 。 在 处 理 非 平 稳 性 及 自 相 似 性 方面 ， 小 波 分 析 
已 被 证 明 是 一 个 非常 有 效 的 工具 。 因 此 ， 在 估算 时 变 Hurst 指数 的 过 程 中 ， 小 波 
变换 将 起 主角 作用 。 对 于 一 个 局 部 自 相 似 过 程 ， 我 们 借助 小 波 变化 建立 了 随机 变 
量 连续 小 波 变换 平方 和 小 波 变换 尺度 之 间 的 一 个 局 部 对 数 线性 关系 ， 然 后 利用 最 
小 二 乘法 估计 出 五 中 。 当 然 在 实际 应 用 中 ， 这 样 的 计算 过 程 是 用 样本 数据 的 离散 
小 波 变换 完成 的 。 


本 章 所 做 的 主要 工作 有 以 下 几 个 方面 : 在 第 二 节 里 ， 在 统计 目 相 似 定义 的 基 
础 上 给 出 了 一 个 局 部 自 相 似 随机 过 程 定义 ， 并 通过 时 变 指 数 的 引入 对 原 有 的 FBM 
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模型 和 ARFIM 模型 加 以 改造 ， 使 之 成 为 具有 局 部 自 相似 性 的 随机 过 程 ， 在 第 三 
节 ， 利 用 Daubechies 小 波 对 局 部 自 相 似 过 程 加 以 变换 处 理 ， 通 过 最 小 二 乘法 建立 
了 一 个 时 变 Hurst 指数 的 估计 式 ， 并 且 证 明了 估计 指数 和 原 指数 是 相 容 性 的 。 在 
随后 的 第 四 节 ， 文 中 给 出 了 获得 时 变 指 数 技术 方面 的 算法 。 最 后 的 第 四 、 五 市 ， 
用 仿真 算 例 对 方法 的 有 效 性 加 以 验证 ， 并 以 上 海 股 市 股票 周 指数 走势 分 析 为 例 进 
行 实证 分 析 。 
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4. 2. 1 一 般 自 相似 随机 过 程 的 讨论 


一 个 统计 自 相 似 过 程 可 大 致 看 成 一 个 和 观测 距离 无 关 且 保持 相同 行为 样本 轨 
道 的 随机 过 程 。 借 助 这 一 过 程 ， 人 们 就 可 以 为 产生 于 许多 科学 领域 中 的 数据 建 模 。 
最 初 自 相似 过 程 是 通过 B.B.Mandelbrot 咎 的 工作 引入 的 。 

设 工 为 随机 过 程 Y() 的 联合 分 布 ， 车 满足 LY(et))= Ze”Y())，c >0， 则 称 
Y(1) 是 一 个 统计 自 相 似 性 的 随机 过 程 。 人 们 在 实际 应 用 中 都 假设 Y() 具 有 二 阶 算 
EL et x (0) = 了 ()-~ 了 Y(t 一 1) 也 是 平稳 的 。 其 中 石 为 自 相 似 参数 ，H e (0,1)。 

假设 了 (1) 为 一 零 均 值 的 随机 过 程 ， 其 协 方差 函数 为 T(s,,s,)】 

D (s,,5,) » CovíY (t ^ s, )Y(t +s.) 


tiselo OR 3E T (s, s, ) 满 足下 式 ， 则 称 F(D) 是 一 个 局 部 自 相似 过 程 ， 
ro)=T00-aGls -s.|" eo), — (4,2,1) 
其 中 g(1}>0， 五 (为 局 部 尺度 指数 函数 。 


由 式 (4,2,1) 式 ， 我 们 可 以 知道 了 (1) 的 数学 期 望 满足 
E{y(t+s,)-Y(e+s, IF f= a(t}s, -s 1 *o(1)j (sil+ls2| っ 0) 


尿 (4.2.1) 式 近 可 知道 , NFRD s, 和 s。, Y(re s) Y (t 4 s, FARE 
因子 js - s," mss. TAER -sl SAQA 1 减 小 到 0 时 ， 随 机 过 程 
7 人) 的 局 部 自 相 关 也 相应 减弱 ， 从 而 产生 出 一 个 逐渐 增 大 偏差 和 粗糙 的 样本 轨道 。 

设 y() 的 增 量 过 程 为 XC)=Y(@)- 了 YLG-1， 该 增 量 过 程 的 协 方差 函数 表达 为 : 


y,(u,,u,)= Covix (t tut +u) 


第 四 章 ” 时 变 分 形 维 数 的 概念 与 计算 m 


假设 增 量 过 程 的 Winger-Ville 4) BAH 7 (A) 
の = ジル クーク "gz 
THA 
A = CO od}, A20, a (4,2,2) 
KH CO) AEA HR. 
当 XAF, 大 (4) 不 依赖 于 !， 从 而 还 原 为 普通 的 谱 密 度 函 数 (A) 


fla)=ca f+} , 420, — een (4,2,3) 


4.2.2 FBM 模型 和 ARFIMA 模型 的 拓展 
1. FBM 模型 的 拓展 : 


我 们 在 下 面 的 一 个 积分 型 随机 过 程 (简称 GFBMO 引入 时 变 的 自 相 似 参数 
H(t)， 可 以 看 成 是 分 数 布朗 运动 (FBM) 的 一 个 拓展 。 


r()- f -O - cuyas) [i-u anu) (4.2.4) 
其 中 実数 7> 0, Blt) od Bus zh. H(n0)e(01). MHRA BR BH 
H(t), GFBM #2 Path ay 35 enc e420, POP) SUUS BLT A. 
FBM 模型 中 的 增 量 过 程 就 是 一 个 分 数 高 斯 噪音 ， 由 于 其 谱 密度 满足 (4,2.3) 式 
要 求 ， 所 以 它 具 有 平稳 性 。 在 通常 情况 下 ， 人 们 总 是 通过 改变 上 自 相 似 参数 鼠 值 的 
方式 构造 出 一 个 改变 了 粗 度 的 样本 轨道 。 相 比 之 下 ，GFBM 模型 的 增 量 过 程 谱 密 
度 函 数 则 满足 (4,2,2} 式 ， 这 表明 增 最 过 程 是 非 平稳 的 。 


2. ARFIMA 模型 的 拓展 


在 分 整 自 回归 滑 动 平均 模型 中 , 我 们 引入 时 变 指 数 右 ((), 对 其 做 一 个 扩展 ( 简 
称 GARFIMA )。 
设 H(t)e (0,1); e(t) Æ Gauss 白 噪音 ，B 为 后 移 算 子 : BX()= BX( -1), 
在 :>0 时 ， 建 立 GARFIMA 模型 如 下 : 
(i - (2) - 8)" ^^ x()- (0 -e(B)e() 
其 中 中 (B)、6(B) 分 别 为 关于 算 子 如 的 PP 阶 和 9 阶 多 项 式 ， 其 特征 根 位 于 单位 
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园 之 外 。 

DP(B)=1-®,(B)- PB’ -----®,B? 
6(B) -1—-0,(B)- 0,B? —...— 0, B* 


另外 一 种 变形 形式 为 : 
o(BXx()-w)- e(B - BY e(t) = (Bn) 
基 中 えり 的 平均 値 。 7( り = (0-8) *"0«(), Hult x() ATE A Eh. “A 


变 分 数 差分 噪声 ”7 人 (导出 的 GARFIMA 过 程 。 
当中 (8)= 9(8) 時 , GARFIMA 模型 可 转化 为 下 面 的 一 个 特殊 形式 .: 
(1- 8) xX()= el) 


考虑 到 等 式 右 端 s(?) 为 高 斯 白 噪音 ， 因 此 这 个 过 程 的 扩展 也 可 看 成 分 整 噪音 


FIN(Fractionally Integrated Noise) 的 一 个 扩展 式 。 

GARFIMA 模型 和 ARFIMA 模型 不 同 ， 它 是 一 个 非 平稳 的 随机 过 程 。 一 个 
ARFIMA 过 程 可 看 成 一 个 具有 全 局 自 相似 性 的 增 量 过 程 目 ， 而 GARFIMA 过 程 则 
是 一 个 具有 局 部 自 相 似 性 的 增 量 过 程 。 在 描述 易 变 形 或 不 规则 性 的 非 平 稳 时 间 序 
列 的 建 模 中 ，GARFIMA 模型 具有 更 广泛 应 用 的 潜力 。 


4.3 利用 小 波 变换 估算 时 变 指数 
4. 3. 1 基于 小 波 变换 的 时 变 指数 估计 


i w(x) Daubechies 母 小 波 ,以 7(a. 忆 代表 局 部 自 相似 过 程 MARTRE 
4a 和 位 置 ! 的 小 波 变换 
WY(a,t)= a^ fy( Y (u)du = a^ jw axes 
由 (4,2,1) 式 及 上 式 可 知 
Ely (at) )= a [fv y eeir eo opua 


= a [f(xy ELV + ax) (e + ay exa» 
= a [fye y OY, (ax, ay )dxdy 
~ a ff ev), (6.0) - atrjax - of ydd 
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zx 4280) (a0) 6e (4,3,1) 
Ext c, =) ff - yd 


Ax y, (a) = logWY (a. t) 


tes fe Ja" oe | 


ya- tegere +c, m 0 om (4,3,2) 
EXC, = E}oghWY(a,t) / E(WY (aty ) 
在 債 祖 小 時 , 由 (4.3.) 和 (4.3.2) 柄 式 可 近 似 地 得 到 一 回 归 模 型 : 
y, (a) (logC,+C,)+2H@) + Mog*+e(a) — eee (4,3,3) 
构造 一 个 小 尺度 序列 ga >a, > >a, a; 227, j 712, ne 
4 x, =loga y, 7 » a, j e 2s ns WF (cy i L2. n] 利用 最小 
二乗 法 得 到 (4,3.3) 式 中 HG eit: 


Alt)= EM - | "m (4,3,4) 


spel, Ye 
4. 3. 2 指数 估计 值 序 () 和 真实 值 五 人 之 间 的 一 致 相 容 性 证 明 


可 以 证 明 上 一 节 的 估计 式 育 ()， 在 n>w 时， 有 EON: 
Var AJ 0 成 立 ， 即 两 者 是 一 致 相 容 性 的 。 
Hi RINER EU) a) 


Y (s, = z) = (log? ph ー n(n * VA K (log? y n(n + aM 一 リノ as 


j= 


> (og 


jal 


= 分 形 理论 中 若干 问题 的 小 波 解 法 
EH (4,3,1) 式 知 : 
loglelwy ay’ je logC, + (2H(t)+ 1)loga + o(1) (a > 0} 
再 结合 (4,3,2) 式 可 推 知 : 
Ey,)= QR (n) 1), +(C, «og? )+ ol) 


将 以 上 各 式 代入 (2.4) 式 中 得 到 ， EAO- | Ese cna 


= H(t)+o(n" m H(t) (n 一 oo] 
其 次 ， 证 明 or 位 全 全 0 
Ti AWM, VarlA())= A eee …(43.5) 


其 中 分 子 部 分 展开 为 
Y, (x, - xs, -Xorb y, )= È, x,x,Covly,,y, )- E>, (x, +x, )covly,.y,) 


+(x)! Y, co) 
设 w AAA M 阶 消失 和 矩 的 Daubechies 小 波 ， 玉 是 一 个 依赖 而 和 a 无 关 
的 常数 。 当 a ,a, 二 0 时 ， 有 下 式 成 并 。 
co mss) = r= fj^ Sy ee 
= [frey OEE + ax) (+ a, yhdrdy 
= (aa, ^ [fye WOT, (21x07 )axdy 


- Jaa, jfr Oe, (0.0) - crai — as" jy 
-[-c 0e | ff la -a an 


= -C(t)a," (かめ a,” Ze + DA yon dudy 


4a, — 0,a, 208^, ->0 时 ， 我 们 将 vu» v, eruditi Mr Taylor 
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展开 ， 并 代入 上 式 得 到 


> + yy, =wlu)t+ yy oA Red Va v Xv, ) {1 o()J 


CovlWY (a, DWY(a,, t)} ~ Ka 05 a, Z ) 
这 个 结果 也 暗示 小 波 变 换 系 数 之 间 的 相关 系数 满足 : 当 被 选择 的 a,,a, 均 很 小 


BA a, <a, RIL 


(ems ] - (4,3,6) 


同时 ， 我 们 知道 ， 对 一 双 零 均值 的 正 态 分 布 随机 变量 X, 了 而 言 ， 有 下 面 的 关 
系 式 成立 : 
lollos(X? Mos?) slo rE om (43.7) 


把 (4,3,6) 和 (4,3,6) 两 式 结合 起 来 得 到 


2(M -H) 


lei fa oes s lot) s. P < (7 ) 


2(M -H) 
于 是 当 a, <a, Ba. a ite, |Cov(rla,).y(a,)) s K GA ) 


由 此 推 知 ; ICov(y,. yj) < d 9 A7 M) 


È, ot Covly,, y j | <K > ijg e-m ) 


< KY ij2 2(j3iXH - M) 


ISj 


=K CM 
i=] j=i 


k 
1992 MEN の (4,3,8) 
im] 


同 理 可 证 
iD, x,Covly,, y 」 ~k’, EYS, Covly,, y, | "m (4,3,9) 


综合 以 上 (4,3,7) 和 (4,3,8) 式 得 到 : 
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Y, -ae -ao 
最 后 ， 将 上 式 代入 Var( 开 () 的 计算 式 (4.3,5) 式 ， 可 以 得 知 


Var(H(t)) ~ k^ — 0,(k > oo) 
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现在 ， 在 离散 、 等 时 空 点 上 观察 某 一 随机 过 程 y()。 

不 妨 将 时 间 点 限制 在 0. 区间 上 ,样本 大 小 为 2 二 =CG-U/a， n=1,2,-°,2°, 
y ak =0 2, 7 204,7 -1) 9 YQ, x27) t8 —A AL) 后 者 基 小 波 変 
换 WY(a,?) 在 a =2-/、t = 天 2 处 的 离散 取 值 ， 采 用 带 有 M YA RAI Daubechies 
紧 支 小 波 基 进 行 小 波 变换 。 

Stepl: 划分 加 ,1 为 2' 个 互 不 相交 且 等 长 的 子 区 间 了,。 

I, -(m-17 m2" is < (J - 1m 212,27 

Step2. 把 应 () 看 成 如 (在 相应 子 区 间 1, 上 的 平均 値 , AO) erg HUE AM 
在 1, 的 中 点 2 (2m 一 忆 。 

定义 如 下 双 变 量 集 : 

(x, ,了 小 = (lose Jiogly sl Ja eJ, O<k<2/-10<j<J- 1 


在 每 一 个 上 上 ， 由 (2.4) 式 估算 H()。 
Step}: 采用 局 部 多 项 式 来 光滑 估算 值 H()， 形 成 一 条 曲线 。 该 曲线 可 看 成 是 
Scop HOER TORT. 


45 仿真 模拟 


下 面 我 们 对 GFBM MGFIN 两 个 模型 利用 上 述 算法 进行 仿真 模拟 。 离 散 小 波 
“As We SEI ELT 8 阶 消失 矩 的 Daubechies 最 小 非 对 称 紧 支 小 波 基 。 选 取 Hurst 指数 
EE H()-3(r-05) +0.1, 1232, J=12, n22^. 


BI 1: HX GFBM 模型 (4,2,4) 式 中 的 B, (t, = 二 六 et) =ー ジ 
n k=) 


i=1,2,…,n; B L) 可 看 成 一 个 布朗 运动 。 对 于 (4,2,4) 式 中 两 个 积分 用 下 边 的 和 式 
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进行 逼近 。 
1 ぐ HUEY HV, 1 ぐ HWY, 
Y (= 一 一 >》 证 一 22 -(- 2 pfu, += (t ~t 2 g(r 
(0- È he -u A Cu p) 1 2-6 7n)" P) 


—— HB o m 
e(u,). e(t JOA Guass 白 品 音 过 程 ，N>0 为 整数 ， 天 >0 为 常数 。 


该 例 GFBM 的 样本 轨道 曲线 见 图 4.1 所 示 。 图 4.2 是 估计 HORSE HO 
线 的 对 比 图 ， 从 图 上 可 以 看 出 万 (信和 估计 羡 ( 人 0 平滑 后 的 曲线 两 着 基本 吻合 。 
例 2, 对 于 GEIN 模型 ， 我 们 采用 下 面 的 一 个 有 限 和 式 逮 近 。 


y i=l 
x)= dab kl) (s UM , 212-4 


Hh et AGRE, Real, QA FATE", 


"P r+ Hl )- 4) 
Å r(e+irlat,)- V) 


部 分 和 式 Y, (r)- VFX, (a) I^ EI E UE OBERE 43 为 
k=} 


和 式 的 一 个 样本 轨道 。 图 4.4 为 估计 AMER HAGR, WALRUS. 


4.6 时 变 Hurst 指数 在 股票 随机 游 动 中 的 应 用 


股票 市 场 是 一 个 自由 度 极 大 的 信息 系统 ， 在 这 一 市 场 中 证 实 一 个 股价 运动 趋 
势 所 需 的 确认 信息 的 时 间 是 不 同 的 ， 对 于 信息 的 不 均等 的 消化 会 导致 一 个 有 偏 随 
机 游 动 。Hurst 和 Mandebrot 俩 人 在 二 十 世纪 七 十 年 代 先 后 对 这 种 现象 作 了 全 面 深 
入 的 研究 ， 并 称 其 为 分 数 布朗 运动 。 这 种 分 形 特 征 结构 的 市 场 实 际 上 会 让 现今 许 
多 诸如 资本 资产 定价 模型 CAPM)、 套 利 定价 模型 APT)、 期 权 定价 模型 等 数 
学 模型 (依赖 正 态 分 布 或 有 限 方 差 ) 屡次 失 录 ， 并 赋予 投资 者 市 场 定性 决策 方面 
循环 、 趋 势 等 概念 及 可 以 度量 的 数量 特征 。 其 中 Hurst 指数 这 一 分 形 统 计量 的 强 
健 性 对 于 时 间 序 列 分 析 分 类 有 着 重要 的 用 途 ， 即 使 是 非 Guass 型 的 随机 分 布 也 是 
可 行 的 。 该 指数 刻 化 了 时 间 序 列 自 相 似 性 及 序列 发 展 的 相关 强度 。 但 是 以 前 人 们 
大 量 的 工作 是 假定 该 参数 是 一 个 常数 。 实 际 上 ， 随 着 市 场 大 系统 非 线 性 演化 ， 过 
程 的 局 部 自 相 似 行 为 会 更 多 地 体现 出 来 ， 这 就 要 求 相 应 的 Hurst 指数 应 该 是 一 时 
变 函 数 H(t)。 与 此 同时 ， 系 统 相 关 增 量 过 程 一 般 也 是 非 平稳 的 。 小 波 分 析 方 法 在 
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此 发 挥 作用 显然 是 理所当然 的 事情 。 

随 着 国家 “ 九 五 ”重大 项 目 《 金 融 数学 、 金 融 工 程 、 金 融 管 理 》 研 究 的 进 一 
步 展 开 ， 小 波 分 析 将 会 在 经 济 、 金 融 领 域 发 挥 重 要 的 作用 。 

本 节选 取 从 1993 年 1 月 1 日 至 2000 年 11 月 31 日 上 海 股票 市 场 周 收盘 指数 
的 所 有 数据 为 样本 ， 按 照 本 文 小 波 处 理 算法 计算 其 时 变 Hurst 指数 ， 并 和 前 人 所 
用 ARFIMA 模型 斗 算出 的 互 ( 常 数 値 ) 相 比较 上 io ， 从 中 可 以 更 清楚 地 看 出 时 
変 Hurst 指数 的 演变 在 投资 策略 方面 具有 非常 重要 的 指导 意义 。 

上 证 周 收盘 指数 共有 398 个 数据 ， 数 据 变化 范围 为 B33.92,2093.2]， 全 距 为 
1759.3 。 数 据 波 动 较 大 ， 旦 现 一 尖峰 分 布 〈 非 高 斯 分 布 )。 | 

当 看 成 常数 时 ， 由 TSP 软件 分 析 ， 得 到 H=0.47， 这 表明 自 1993 年 以 来 周 指 
数 之 间 波 动 没有 明显 的 长 期 依赖 性 。 利 用 我 们 的 算法 及 三 次 样 条 曲线 拟 合 ， 从 而 
得 到 时 变 Hurst 指数 的 曲线 ， 可 以 明显 地 看 到 ， 在 股市 不 同 波段 中 周 指数 依赖 性 
都 有 不 同 的 变化 ， 这 对 实际 的 投资 决策 有 着 潜在 的 影响 是 不 言 而 渝 的 。 

图 4.5 为 上 证 周 指数 折线 图 ， 图 4.6 则 显示 了 时 变 Hurst 指数 变化 曲线 及 传统 
方法 计算 的 指数 值 图 线 。 
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在 这 一 章 ， 我 们 提出 了 时 变 维 数 D() 的 概念 ， 通 过 分 析 时 变 Hurst 指数 万 (のり 
来 实现 这 一 目的 (DQ)= 2 H(f))。 在 定义 了 统计 局 部 自 相似 过 程 之 后 ， 将 时 变 指 
数 玉 从) 分 别 引 入 分 数 布朗 运动 (FBM) 模 型 和 分 整 自 回 归 滑 动 平 均 (ARFIM) 模 型 之 
中 ， 对 它们 各 自 进 行 了 拓展 ， 成 为 典型 的 具有 统计 局 部 自 相似 性 的 随机 过 程 。 

由 于 拓展 后 的 两 个 模型 均 具 有 非 平 稳 的 性 质 ， 文 中 借助 Daubechies 紫 文 小波 
变换 ， 给 出 了 时 变 指 数 估 值 良 (1) 的 小 波 算法 ， 并 证 明了 估 值 8() 和 真实 值 (0) 是 
一 致 相 容 性 的 。 

最 后 ， 对 于 两 个 模型 分 别 进行 了 仿真 模拟 ， 结 果 比 较 理想 。 另 外 ， 通 过 上 海 
股价 周 收盘 指数 的 时 变 指 数 产 (f) 的 计算 分 析 ， 本 文 认为 ， 时 变 指数 (0) 的 引入 比 
原来 常数 指数 互 对 于 股票 投资 决策 来 说 具有 更 好 的 指导 作用 。 
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分 形 最 突出 的 特征 是 其 在 不 同 尺度 下 所 表现 的 自 相似 性 。 统 计 自 相似 随机 过 
程 则 是 在 概率 分 布 (g) 意 义 下 满足 自 相 似 性 ， 
faa" sf (at) 


在 这 个 定义 下 ， 我 们 可 以 看 到 这 个 过 程 在 均值 、 自 相关 函数 及 功率 谱 方面 均 
呈现 出 自 相似 性 ， 所 以 也 就 呈现 出 了 分 形 的 分 维特 征 。 
所 谓 /过程 就 是 一 类 功率 谱 具 有 如 下 形式 的 随机 过 程 ; 


人 2 市 | 


FT 
S loj=— 


现实 生活 中 ， 诸 如 真空 及 电子 元 器 件 中 的 噪音 、 股 票 指数 的 变化 、 生 物 心率 
跳动 、 网 络 流量 等 大 量 的 自然 现象 都 具有 这 类 过 程 的 特征 。 

对 这 种 Vo 辻 程 的 刻 画 , Mandelbrot 等 人 从 时 域 上 提出 分 数 布朗 运动 模型 
(FBM) 和 分 数 高 斯 噪声 模型 (FGN) 均 基 本 上 可 以 解释 这 一 过 程 ， 但 总 的 说 来 ， 出 
T Y 过 程 的 非 平稳 性 ， 其 处 理 方法 较为 复杂 ， 对 谱 参 数 y 的 选择 也 不 近 人 意 。 
与 此 同时 ， 过 程 的 非 平稳 性 也 为 小 波 分 析 方法 进入 这 一 领域 提供 了 合适 的 场所 和 
机 会 。1992 年 G.Wornelko] 给 出 了 刻画 多 过程 的 小 波 模 型 摆脱 了 FBM 模型 的 


局 限 性 。 有 关 V 过 程 谱 参 数 估计 问题 ， 随 后 他 提出 了 用 极 大 似 然 估 计 方法 估算 


”参数 的 途径 Du23 。 
作 极 大 似 然 估计 要 从 观测 值 的 似 然 函 数 ， 也 就 是 从 条 件 先 验 概率 入 于 。 
BERMA- A x() RRS w(x) RR: rn) xt x(r)-w(). EKER 
度 六 下 的 小 波 变换 后 六 ” = ai? tw", BA Varr. と "| = o? =0 p" +0}, B-2. 
G.Wornel1992 EE XP Hir" ET Gauss 分 布 ， 得 出 似 然 函 数 为 : 


de 
pA)= n e 


mnm mt foro? Jro} 


RT BAL A BM A 


Ô; i Fs») Sy WE RE m 下 的 样本 方差 。 
ke N(m) 


xt (5,11) 式 分 別 対 ,o?,o25R S, HARM SION E Af 68:801 3EÉX 


性 方程 组 : 


沪 方 程 组 的 求解 可 用 文献 2 中 提出 的 数值 方法 加 以 解决 。 
从 以 上 过 程 可 以 看 出 ， 这 种 方法 要 求人 们 在 使 用 该 法 之 前 要 知道 随机 过 程 的 
概率 分 布 函数 ， 应 用 起 来 带 有 许多 分 布 假设 的 成 分 ， 这 是 其 不 足 之 处 。 


通过 对 V, 过程 的 小 波 变 换 特 征 的 研究 ， 本 文 欲 借助 小 波 变换 对 统计 自 相似 


过 程 所 起 的 白化 作用 及 最 小 二 乘法 ， 建 立 起 在 无 须知 道 样本 概率 分 布 情况 下 的 谱 
参数 估计 方法 。 

本 章 所 做 的 具体 工作 是 围绕 以 下 几 个 方面 展开 的 ; 

在 第 二 节 里 ， 文 中 探讨 了 V 过 程 的 频 域 表 征 及 其 非 平稳 性 ， 随 后 ， 在 第 三 


节 着 重 就 有 关 V 过 程 小 波 变换 系数 之 间 的 尺度 相关 性 及 小 波 函 数 消失 矩 对 系数 


相关 性 衰减 的 影响 进行 了 详细 的 讨论 ; 基于 最 小 二 乘法 的 谱 参 数 的 估算 安排 在 第 
四 节 进 行 ， 文 章 的 最 后 用 算 例 加 以 模拟 验证 。 


S2 YY 过 程 的 频 域 表征 及 其 非 平稳 性 


= 分 形 理论 中 者 干 问题 的 小 波 解法 


设 一 个 统计 自 相 似 过 程 x(/) ( 功 率 有限 ) 的 自 相关 函数 为 Rs(r) DRWA 
S,(w)， 其 相互 之 间 的 关系 是 : 


+T 
R,(c)= lim ーー O+, Rr) = an, (ar) 
zT 


S,() = [R, (ce de =a 9"? fR, (ry "ar (r =ar) 


ー の -(2H “Ne (2) 
a 


= a?" S (aa) 
由 此 可 以 推 知 ， so), prar Cro lona. 


4 y=2H+1, W sos Hoy y 为 自 相似 参数 。 


上 式 表 明 x(1) 的 功率 谱 具 有 分 维 谱 特 征 。 
S.(@) 的 非 平稳 性 体现 在 以 下 两 个 方面 一 方面 ， 当 y>1 时 ，w 愈 小 ， 则 


S_(@) 就 会 愈 大 ， 说 明 低频 成 分 丰富 ， 但 当中 一 0 时 ，S.(w) 一 ， 从 而 造成 所 
iBh) "£L AGE", 男 一 方 面 , 当 ヶ <1 時 , の 候 大 , US (o) EC. Mo — o 
時 BAAS (0) > o ， 又 造成 所 谓 的 “紫外 灾难 ”。 为 了 避免 这 种 两 极 现象 的 发 
生 ， 有 必要 在 频 域 上 采取 一 些 限 制 条 件 ， 对 过 程 加 以 定义 。 

X521 一 个 广义 自 相似 的 随机 过 程 称 为 一 个 VL 过程 是 指 ; 


—m | jl e,«o«o, 
设 有 理想 带 通 滤波 器 if(o)- ite: 

如果 存在 の 。 和 の (0 < の 。 <@w <o)， 使 得 当 x( 作 为 输入 过 程 通过 一 个 具有 
频率 响应 的 理想 带 通 滤波 器 访 (w) 后 ， 其 输出 y(f) 是 具有 有 限 方差 的 广义 平稳 随机 


可 以 证 明 ， Am 


" 形式 的 VER x0) & E35 Hio) MEER, HEPES UY 
e 


稳 的 随机 过 程 yf) ， 其 功率 谱 S,(@) 满 足 : 
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< wl < oy 


9 
3, (a) = lol" 
0 AUC 


5.3 V. 过程 的 小 波 变 换 在 相关 结构 上 的 特点 


Bal ERREA ;的 上 过 程 W(t) A RO REND RB xl) 
在 二 进 制 小 波 变换 下 的 系数 为 dl”'， 即 就 是 


ftv Oat =0, (o0«&k«R), | d? -(x()w,,)-2* | ylane- nyit 


设 采 样 数 率 为 T ， 对 某 一 信号 x(f) 的 采样 序列 a? = x(nT,)， 可 看 作 是 尺度 
カ =0 的 対 x( り 的 近似 。 当 > 0 时 ， A RR df") 之 间 有 下 面 的 递 推 关系 式 ; 
a T Dhow. T 


ae zi > Bsp" 
ERP h 和 gg 分别 为 尺度 函数 p(x) 和 小 波 函 数 y(x) 对 应 的 滤波 系数 。 
由 Mallat 重 构 算法 得 到 : 
a, = S ha + Epu de 
下 面 探讨 任意 两 个 小 波 变换 系数 do 和 df] 之 间 的 相关 性 : 
Ela! mgl |= N Ü Ely,,, Oxel y, , (r)htar 
=[ Ev. Rt ypy € utar 
= [Vn OR ry OF 
利用 Parcival 定理 可 知 上 式 可 化 为 


glgege) に uu era pr ab peep 


27 
hem] noon 
=e [zw mo 外 p "m COT Pde --(5.3.1) 


20 lg 


由 上 式 可 得 到 以 下 结论 : 


分 形 理 疹 中 若干 向 是 的 小波 解法 
|. Da RA” OLA EB, DA Varda o2" 。 


由 于 小 波 函 数 y() 满 足 [yur = 0 ， 因 此 小 波 系数 dim) 的 均值 为 零 。 
TE (5.3.1) P m e m,nsn ， 于 是 方差 为 : 


Ela P = varla] -= Z- is S ("oJ do 
ERU XM 0 iussi (5.3.2) 
ES 
这 里 ， c! = う ァ “ip " (o) do 
利用 这 一 


结论 将 (5.3.U 式 进行 归 一 化 ， 归 一 化 的 相关 系数 为 : 
Ea ae) 
Var(d Wara) | 


tytn) wong 


» : di 27" 9 see je. 
y T y Ed 


: @ 
の の 


PELA 108); CE 
2 ド SAC LB, cos|n27" =n" o to 
び を 0 a’ 


p(n,n';m,m )- 


ト ュ 


SS a P ene ee 
在 (5.3.3) 式 中 , mom, 得 到 


o: の -y) *o 


pn n ;m)- pu E piai | cos|2 m(n n oho 


=" eer cos|(n — n Ho 


上 式 表明 : 在 固定 尺度 下 ， 小 波 系数 d 人 中 是 广义 平稳 的 随机 序列 。 
3. 在 不 同 尺度 下 同一 时 刻 的 小 波 变换 ， 归 一 化 相关 系数 为 (mm ) 的 函数 在 
(5.3.3) 式 中 ， 令 12-” =n2-”， 这 样 就 得 到 


Mme m Yi- ) + -m + nm 
pln,n ;m,m )= a ーーー に = < dw 
0 
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2 Men! \i-y) toy -(m-m } * 
o? 2 pre @ (0) 7 


o T : c 


该 式 表明 : onn mm ) 与 采样 时 刻 的 值 没有 关系 。 
4. 小波 函数 A 際 消失 逢 対 小波 系 数 相 共 性 的 影 咽 。 
(D A Var(at oy A BEL, RUTMER aE 


ーーー 
— 


正明 : HH, w(e)-o'w(»), (¥,(0) 0) 


o? +0 ny 
g!- 2c Sap ~ by, (oy do 


0. 12A- 2 
= fel") do 
MM 


上 式 中 内 有 当 7 <2R 时 ，o? 才 是 有 限 的 。 
(2) 在 同一 尺度 下 ， 归 一 系数 p 值 随 着 (nn ) 增 大 而 衰减 ， 这 种 衰减 量 
Bh (n - pi a 


m 2 
pln, n ;m,m )- : - [EEL ofn -r pho 
0 


g? 1 2R-y 

gp lla lb, (o) cos|(n - n ho 
AERA, or nm) BS (n-n CK RUN, 3E 

REA (n-njo Eee. MAAS, PR BY RR R 

大 ， 相 关系 数 p 衰减 就 越 快 。 这 一 结论 给 我 们 最 大 的 启示 是 ， 可 以 

近似 地 将 小 波 变换 系数 4 人 "看 成 是 一 个 白 噪声 序列 。 


5.4 ”基于 小 波 系 数 方差 结构 的 y 参数 估算 


过 程 的 一 个 样本 函数 x(1) 进 行 采样 ， 利 用 离散 小 


对 于 具有 功率 谱 ， 
e 


波 变 换 处 理 。 不 妨 设 选用 的 膨胀 尺度 为 m ， 对 应 的 小 波 系 数 个 数 为 N(m) 。 
由 于 小 波 系 数 之 间 的 相关 性 较 小 ， 可 以 忽略 不 计 。 对 于 尺度 m ， 小 波 系 数 的 
方差 Jarle 四 | 可 以 估计 为 : 
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XH (5.3.2) 式 知 : rlg の |= g?” 

建立 回归 方程 : G? =0?2" 

PR IZ. BOSE Ae 7S BUF LII E PE Ar Fe 
fn = 8- my 


ws f= log" ， g=log? 


b difti) 


n=l 


TÉ g Aly B — Sei EE Prahe y. 
即 P(8,7)= Min P(g.y) 


由 多 元 函数 极 值 定 理 知 ， 少 和 & 是 下 列 方程 的 解 : 


as 7 in B+ mr)-o 


m= 


o, Dn, -gt my)=0 
功率 谱 参 数 的 估计 值 为 : 


73 mf, + 6(M + D f, 


m=] m=} 


oe 


55 算 例 仿真 


对 于 y =1.67 的 过程 ( 见 图 5.1)， 采 集 样本 个 数 N = 2" ， 选 用 具有 4 阶 


HABI] Daubechies 小 波 基 估 算 7， 仿 真 结果 见 下 表 5.1 及 图 5.2 所 示 。 
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# 5.1 
估 值 7》 绝对 误差 相对 误差 
DT E 1. 65 [ 20.02... | L% | 


从 结果 来 看 ， 本 文中 的 方法 是 有 效 的 。 


5.6 ”本 章 小 结 


在 本 章 中 ， 我 们 首先 在 频 域 中 描述 了 V 过程 的 特征 ， 并 着 重 探讨 了 Vii 
程 经 过 小 波 变换 后 小 波 系数 在 相关 结构 上 的 特点 ， 有 关 结 论 主要 有 
(1) 在 小 波 的 同一 尺度 m 变换 下 ， 归 一 化 相关 系数 是 (nn ) 的 函数 ， 表 明 
小 波 系数 4 中 是 一 个 广义 平稳 的 随机 序列 。 
(2) 在 不同 変換 尺度 下 的 同一 時 刻 Caa" =n2-* )， 归 一 化 相关 系数 为 
(m-m ) 的 函数 ， 与 采样 时 刻 的 值 没 有 关系 ， 
(3) 小 波 的 消失 矩 越 大 ， 小 波 系数 的 相关 性 衰减 就 越 快 。 
关于 V 过 程 谱 参 数 y 的 估计 问题 基于) 过程 的 小 波 模型 ，G.Wornell 提 
出 了 用 极 大 似 然 方 法 估算 y 参数 的 途径 ， 但 存在 着 须 预先 知道 过 程 的 分 布 函 数 的 
不 足 之 处 。 从 上 述 的 相关 性 分 析 可 以 看 出 ， 小 波 变 换 对 于 自 相似 过 程 具有 白化 作 
用 ， 以 此 为 出 发 点 ， 本 章 提出 了 无 须知 道 分 布 的 情况 下 谱 参数 y 的 一 种 估计 方法 。 
最 后 ， 利 用 一 个 仿真 结果 对 该 方法 加 以 验证 ， 说 明 此 法 是 有 效 的 。 
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图 5.2 估算 值 y = 1.65 
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第 六 章 多 重 分 形 分 维 谱 的 奇异 性 及 其 小 波 分 析 
6.1 引言 


在 研究 复杂 分 形 《〈 如 混沌 现象 等 ) 的 过 程 中 ， 人 人 们 发 现 仅 用 一 个 分 形 维 数 来 描 
述 非 线性 动力 学 演化 过 程 中 产生 的 形体 显然 是 不 够 的 。 不 同 的 局 域 条 件 或 由 参量 涨 
落 的 波动 均 会 造成 形体 形态 方面 的 差异 。 为 了 弄 清楚 在 分 形 形 成 过 程 中 局 部 条 件 的 
作用 及 分 形 人 参量 在 空间 的 分 布 问 题 ， 和 人们 提出 了 多 重 分 形 (multifractaj) 的 概念 (也 出 
分 形 测度 ) 。 早 在 1974 SE Mandelbrot 就 首先 指出 ， 在 流体 发 生 湛 流 时 出 现 了 多 重 标 
度 特性 ，“ 自 然 界 的 分 形 几何 学 ”一 书 中 ， 他 用 “ 非 间 际 的 ”Qnon-lacunary) 以 此 来 
描述 这 一 现象 种。 随后 ，Grassberg 等 人 在 1983 年 研究 奇异 吸引 子 时 发 展 了 这 一 概 
念 ， 并 提出 现今 使 用 的 一 些 数学 公式 。Halsey A M WE 1986 年 应 用 这 一 理论 研究 
了 有 限 扩 散 凝 聚 (DLA) 集 团 的 分 形 生 长 问题 。 我 国学 者 段 钻 荣 和 颜 新 祥 1995 年 5] 06 
尝试 将 多 重 分 形 的 多 标 度 性 与 小 波 变 换 良 好 的 局 部 分 析 性 质 相 结合 ， 在 预测 油气 田 
储量 及 勘探 井 位 的 应 用 方面 获得 了 令 人 满意 的 成 功 。 近 来 ，Rudolif.H. 等 人 提出 了 一 
个 多 重 分 形 的 小 波 模型 来 解决 网 络 流量 问题 po 。 总 之 ， 多 重 分 形 在 解决 实际 问题 中 
发 挥 着 重要 的 作用 ， 并 且 小 波 工具 也 已 经 开始 进入 多 分 形 的 研究 领域 。 

概括 起 来 说 ， 多 重 分 形 外 观 上 看 是 不 同 无 标 度 区 的 集合 ， 它 将 多 个 维 数 汇 成 一 
个 系统 一 一 分 形 谱 函 数 ， 主 要 研究 分 形 测度 《物理 量 ) 在 文集 上 的 分 布 情况 。 县 有 
紧 支 集 的 小 波 在 空间 是 局 部 支撑 的 ， 因 此 可 以 利用 小 波 分 析 解 决 多 重 分 形 谱 的 局 部 
细 化 和 奇异 性 变化 问题 。 在 这 方面 ，Mallat 和 Hwang 在 1992 年 利用 小 波 变换 模 极 大 
值 法 开创 了 研究 信号 〈 包 括 分 形 信号 在 内 ) 的 局 部 奇异 性 行为 的 先河 5 ， 为 以 后 人 
们 专门 研究 多 分 形 谱 的 奇异 性 提供 了 理论 研究 的 线索 和 思路 。 


在 前 人 已 有 多 重 分 形 理论 的 基础 上 ， 本 章 探讨 了 如 何 应 用 小 波 变换 处 理 有 关 多 
分 形 谱 的 奇异 性 问题 ， 其 主要 工作 和 创新 之 处 包括 在 以 下 几 个 方面 : 

首先 在 第 二 节 介绍 了 目前 描述 多 分 形 的 两 套 语言 一 -wx- /oj 语言 和 4- 疡 语言 
及 两 者 的 等 价 性 。 在 第 三 节 中 ， 在 Malat 和 Hwang 用 小 波 模 极 大 值 对 分 形 信号 进 
行 奇 异性 分 析 的 基础 上 上， 以 高 斯 函数 作为 母 小 波 ， 提 出 了 一 种 跟踪 每 一 条 小 波 模 极 
大 值 线 的 算法 ， 从 而 避免 了 大 量 的 数值 运算 。 为 了 更 进一步 解决 奇异 点 附近 可 能 存 
在 的 一 些 拓 扑 问题 ， 文 中 第 四 节 采 用 了 两 个 多 项 式 函 数 作为 芽 函 数 探讨 了 小 波 变换 
极 大 值 线 的 拓扑 结构 。 本 章 第 五 节 提 出 了 基于 小 波 变 换 模 极 大 值 方法 计算 多 重 分 形 
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谱 的 一 种 方法 。 本 章 的 最 后 采用 了 两 个 范例 ， 其 一 用 于 说 明 分 形 测度 奇异 性 和 普通 
信号 奇异 性 的 不 同 之 处 ， 这 是 探讨 跟踪 问题 的 必要 所 在 ， 其 二 用 于 解释 多 分 形 谱 的 
计算 。 


6.2 多 重 分 形 的 两 类 描述 语言 
6. 2. 1 a- ro) 语言 


议 x 为 拓扑 维 数 为 4 的 任 一 支 集 ，4w 是 zx 上 的 一 个 测度 。 对 于 {x, 2) a 4 
式 划分 ，@ 是 和 划分 有 关 的 参数 。 第 n 步 划分 后 具有 相同 测度 jy 的 单元 构成 的 支 集 
A x, (a). 

T lim x, (a) x, 是 一 个 分 形 集 ， 则 称 x, A (x, 4) 的 一 个 分 形 支 集 (如 Cantor 集 
中 的 三 分 Cantor 集 情 形 ) 。 于 是 一 个 多 重 分 形 可 表示 为 具有 不 同 维 数 的 分 形 子 集 x, 
的 并 集 。 

右 单 元 测度 u 与 単元 的 尺度 之 同 満 足 零 律 共 系 , / -g“ ， 则 称 w 为 Elder 指 
数 。 因 为 a 控制 着 测度 jy 的 奇异 性 ，a 也 叫 奇异 性 指数 。 

XFER x, 的 任意 6 - Em (U,], 、， 当 0<diamU, < 65 时， 引入 以 下 记号 : 

H; (x,u, )= int È lant, Y や の Qu. 
类 似 Hausdorff 測 度 的 定 叉 , EX x, 的 7 维 测度 为 : 
H' (Xas )= lim H (xq, Ma) 
若 存在 临界 指数 f(a), [518 
0 r> fla) 


oa oo r< f(a) 
ARES ”r= f(a) 


则 称 /(a) 为 多 重 分 形 在 Hausdorff 测度 下 的 奇异 性 谱 ， 即 就 是 由 不 同 w 组 成 的 序 
列 构成 的 谱 。 
HERE MAT, f(@) 就 是 分 形 支 集 x, 的 Hausdorff 维 数 。 


f(a) = inf I lim int 2 \diamU, y = | 
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= sup lim inf ゞ (diamU, Y = «| 
» ist 


4 ERAZ U) ,是 尺度 为 < 的 盒子 时 ， 若 [a,a+ da| 内 测度 为 u 的 盒子 数 为 


N(a) , WA FARMA: 
Hi(x, Ho )= N(a)e’ 


b ， Hia H HARE. Hert f(a) = mn 


下 面 从 信息 论 的 角度 给 出 另 一 套 描 述 多 重 分 形 的 参量 。 
设 测度 支 集 x 经 过 迭代 或 划分 后 的 单元 为 {4,}， 概 率 测度 为 x 。 第 1 单元 的 概率 


定义 为 P: 
P = | aulx) 


当 ア > アッ Hq»»18j, 有 P* om M. ESL u Wg BEA: 
N(q,«)- 2 P(e) 


Mi(q)= > P (e) e” = N(qe}e" 


M’ (q)= lim M; (g) 


车 存在 依赖 于 P 的 g 阶 和 矩 所 选择 的 临界 指数 r(qg)， 使 得 
0 r > t(q) 


M'(q)=4 œ r < tla) 
有限 正 数 ァ =z(2) 


別称 zr(g) AAEE. H u Hg NIE, 定义 广义 Reny 的 g 次 信息 维 如 下 式 
l InN(g,e) 
pum ye (g #1) 


D, = £-»0 
| nz(9/ (q 21) 
lim 2,7 Ing 


6. 2. 3 a-f(@)Batg-D, IBEBISEDETE 
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D, 和 r(gj 有 如 下 关系 ; 


r(a) z1 
"E | ー mn q = LBT 


若 r(gj 和 oj) 均 可 微 时 ， 可 实施 如 下 的 Legender 变换 ， 
Fle(a)] - ao(a)^ rtg) 


aq) - -7 0) 


p, = o (a)- f 2 A 


a(q)= dl(q - LA 


7@)= aa -|a(a) - 110,6; 
从 上 述 推 导 的 过 程 可 见 ， 在 已 知 g - D, 的 情况 下 ， 另 一 套 参 量 g f(a) 就 可 由 上 
面 两 个 式 子 求 出 。 


于 是 得 到 


63 多 重 分 形 奇异 谱 结 构 的 小 波 分 析 
6. 3. 1 前 人 有 关 奇 异性 的 局 部 刻画 的 工作 


一 个 函数 (分布) 在 某 一 点 x 处 的 奇异 性 可 以 用 其 Holder 指数 h(x ) 来 描述 。 
Ax。) 越 大 , 7(x) 在 x。 邊 就 越 正規 。 若 ヵ < h(xo)<n+1， 那 么 我 们 就 称 f(x) 是 n 次 
(而 非 4+1 次 〉 可 微 。 传 统 方法 是 分 析 f(x) 的 Fourier 変換 (o) 34 c っ o 時 的 河 近 
衰减 情况 。 由 于 Fourier 变换 缺乏 空间 局 部 性 ， 此 时 这 种 渐 近 衰减 只 反映 了 有 关 区 间 
上 发 生 奇 异性 的 整体 “宏观 ”信息 ， 而 难以 确定 其 奇 点 的 空间 分 布 。 
相 比 之 下 ， 小 波 分 析 则 可 提供 有 关 奇 异性 的 局 部 信息 。 
(1) Argoul.F.,Aneodo.A. 和 Bacry 等 人 在 90 年 前 后 用 小 波 方 法 开展 对 分 形 测度 
性 质 进行 分 析 ， 发 现 分 形 测度 的 小 波 变换 与 局 部 尺度 指数 之 间 存 在 有 下 面 的 一 个 帮 
律 关 系 : 
假设 以 x 为 中 心 、& 为 长 度 的 区 间 7(x,,e) 的 分 形 测 度 为 x ， 且 满足 表达 式 : 


ul (x, ,Ae)| = A" urs. e) 
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其 小 波 变换 定义 为 : 
W, (a,b)= a p(t EAS 


其 中 a AMEBIS ve Abie RB y E E eR 
则 有 下 述 结论 成 立 : 


W, (4a, x, + Ab) = NE ^^ lu) 
f の M2 ue ex) & xx, 
"(uy M 3 Pte) 


~are (a,x + の (2 っ 07) 
这 里 alx ) 是 分 形 測度 (>) 在 x, 邊 的 局部 尺度 指数 ( Holder exponent ) 

(2) 为 了 分 析 信 号 局 部 的 奇异 行为 ，Mallat 和 Hwang 引入 小 波 变换 模 极 大 值 
分 解 方法 ， 此 法 可 用 来 定量 描述 分 形 信号 的 奇异 性 四 ] 。 

为 了 刻画 在 某 点 〈 而 不 是 某 区 间 ) 的 奇异 性 ， 设 wy(x) 是 有 紧 支 撑 的 的 n 次 连续 
AY GANT RPA, AAG nA, 有 如 下 定理 : 

定理 : Snen, a(x, ) Av Tn 的 非 整 数值 ， 分 形 测度 函数 u(x)e LR). $ 
u(x) ZE 5 x, 是 Lipschitzg(x。 ) 的 。 史 存在 一 常 数 ん > 0， 使 得 对 于 xo 邻 域内 的 所 有 点 
和 任意 尺度 a ， 分 形 测度 的 小 波 变换 满足 

W,u(x) < Ka?) 


或 存在 常数 K > 0 AMF x, 邻 域 内 的 所 有 点 和 任意 尺度 a ， 分 形 测度 的 小 


波 变换 满足 
a(x) 
ZU e+ 区 -得 | 


logi — xol 


但 在 实际 计算 中 应 用 上 述 定理 是 相当 困难 的 ， 因 为 需 在 U(x E) HI Zt (a,x) Z 
同上 捜索 u(x) 的 衰减 性 ， 除 非 所 有 奇异 性 有 相同 的 指数 a(x。), TU, RFEA 
量 的 计算 量 。 

XIF VxeU(Gs c). A Wulay.x) xWu(ay x). WE (aj, x, ) 为 小 波 变换 横 极 
大 值 点 ; RETN (a,b) a > 0,5 > 0 中 的 一 条 曲线 工 为 小 波 变换 模 极 大 值 线 ， 若 
vc.5)sT， 则 (@.) 为 在 尺度 wa 下 的 模 极 大 值 点 。 
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利用 这 种 小 波 变 换 模 极 大 值 方 法 可 以 较 精 确 地 检测 信号 奇 寞 点 的 位 置 ， 但 对 于 
(E) 分 形 来 说 ， 由 于 奇异 点 是 非 孤 立 的 ， 这 就 很 自然 地 涉及 到 跟踪 问题 ， 这 是 一 
个 空白 点 。 

总 的 说 来 ， 上 述 方 法 的 不 足 之 处 可 归结 为 下 面 丙 点: 其 一 ， 在 处 理 分 形 信号 非 
拆 立 奇 寞 点 时 没有 跟踪 能 力 ， 这 就 意味 着 分 形 在 细节 上 有 部 分 信息 丢失 ， 分 形 谱 的 
计算 就 会 失真 ， 其 二 ， 局 部 极 大 值 是 在 离散 尺度 上 并 且 在 局 部 菜 一 尺度 级 水 平 上 形 
成 极 大 值 线 ， 在 离散 点 (a,, ,6b,,) 上 计算 要 涉及 大 量 的 数值 运算 ， 此 外 也 没有 考虑 极 大 
值 点 附近 的 拓扑 问题 。 

下 面 ， 在 微分 方程 数值 解 的 基础 上 ， 我 们 提出 一 种 寻求 最 优 尺度 (a,5) 来 跟踪 每 
一 条 极 大 值 线 的 算法 。 


6. 3. 2 高 斯 函数 的 一 些 性 质 


为 了 减少 计算 量 ， 了 解 测度 函数 Ap(xz) 中 具有 Holder 指数 hx)<n 的 所 有 奇异 
性 ， 需 要 一 个 具有 nn 阶 消失 和 矩 的 小 波 母 函数 。 在 这 方面 我 们 选择 高 斯 函数 作为 母 小 
波 ， 它 具有 民 好 的 局 域 性 且 容 易 满 足 消失 矩 的 要 求 。 

XE X. n BT B E ER CJ 


v, (x)= y, ax): Volx)=e 


d 

ax 

n>0 — w,()- yw, (x) 
dx" 


n< 0 时 W n (x)= fimi (£e a C. 


xo € 
2 


*C,-20. nz-1, & C(x)= (x)= fe 
于 是 Vne Z A FIBBJXCRSUXV: 


(n+ I, (x) + xv, x) WY (x) = 0 | (6.3.1) 


d$ 


由 此 式 可 推 知 
w(x)= P(x)yo(x)+ 0,(x)C(x)}, (ne2) 


其 中 P(x) 和 の O (x) 均 是 多 项 式 ， 有 下 面 的 递 推 规律 成 立 : 


(n+ 1)P (x) + xP... (x) + Pa(x)=0 
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(n a 1)0, (x)+ xO, 41 (x)+ n+2 (x) -0 


6. 3. 3 跟踪 高 斯 小 波 极 大 值 线 的 参数 法 


"— id NIAE TR MA 
W" (2, の ) = ; J^ (x, z=" Jas (a,beR',a>0) 


d 
BS ga ep 有 2-2 
y2*- b ， 则 "ria p 


gi d* d' : 
id Wn m r a E u) 


由 (6.3.1) 式 可以 推 知 : 
Sal) = Ca vu) 


d' 


i E 
(dw,0) D 2, isa) 


AZ n 2i)- (& - 1)] O<i<sk#0 
其 中 A= 
0 其 它 


= B= -i+1), (k>0) 
于 是 可 得 到 小 波 变 换 的 导数 运算 表达 式 


W' uH ase C Ya y な ia. "s ly awe uu 


1-0 
"Ak =1 时 ， 我 们 得 到 一 个 微分 方程 
aW? u-W' us Wp = リ 0, 
a 


对 于 给 定 的 尺度 a。， 在 b 处 的 局 部 极 大 值 由 下 边 两 式 确定 


Oe fe 
W, u(a,,b) 。 -0, ab? "4 CN b) " «0 


这 些 极 大 值 点 的 连 线 叫 极 大 值 线 。 


pisa 


te 分 形 理论 中 者 干 问题 的 小 波 解 法 


前 人 的 工作 都 是 在 离散 尺度 a MAE b, 情况 下 提取 局 部 极 大 值 ， 然 后 把 这 些 极 
大 值 点 连接 起 来 。 现 在 假定 小 波 变换 是 实数 连续 的 ， 我 们 相信 ， 在 一 个 最 好 尺度 上 
所 有 极 大 值 点 都 可 以 确定 下 来 ， 而 每 一 个 极 大 值 点 均 可 看 成 一 条 极 大 值 线 的 起 始 
扎 。 为 此 将 a、4。 分 别 看 成 是 参数 1 的 函数 alt)、6b(1)， 以 (alt)5() 来 表达 一 条 极 大 
值 线 ， 这 样 沿 该 条 极 大 值 线 极 大 值 点 的 运动 由 下 面 的 一 个 偏 微分 方程 决定 : 


d da db 
AW, u)=Wrus+wrn® =0 
2 "= Uy eA 
da db 
— -—CW'u, 2 EW’ 

等 价 于 P phe ba H 


ZEW! =O RUE F, 可 用 Runge Kutta 方法 求解 上 述 方程 。 
6.4 多 分 形 测度 小 波 变 换 极 大 值 线 的 拓扑 分 析 


当 小 波 变 换 极 大 值 点 附近 比较 平滑 或 着 在 分 贫 点 的 邻 域内 ， 上 述 分 析 方法 的 效 
呆 就 不 太 明 显 ， 因 此 有 必要 探讨 这 些 极 大 值 线 的 拓扑 情况 。 
在 局 部 静态 分 岔 问题 中 ， 关 于 奇异 点 附近 的 函数 性 态 ， 通 常 选用 一 些 简单 的 多 


項 式 画数 GS Hist CGF) sex , 
fW, ula, b) ERW u(a,,b, ) = 0 的 邻 域 内 展开 为 Taylor 级 数 : 


Wy u(a,b) た の De Jaa ADW” gn „dla ba) = T(Aa, Ab), …… (6.4.1) 


i=0 6 


其 中 Aa=a-a,, Ab=b-b, 
iW u Elab Ab 1-18T38 x, 即 就 是 


0 
W^ ulag b, ) = | 0 


l<m<il 
m={ 
这 里 ， 取 两 个 简单 情况 加 以 说 明 : 当 7= 20[. (a,b, ) 就 是 极 大 值 线 上 的 一 正则 
点 ， 而 当 1 =3 时 就 会 有 一 个 分 岔 点 产生 。 
令 = 各 ，1= 望 ， 由 (63.2) 式 和 (6.4.1) 式 可 以 推出 
a a 


I n+l ご I 7 H-(11),2:- 01 
T= W t 
p a ) d "e; E |, 
Ig 3 T p^ ~e} ルポ / y 


a 7=0 
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251 为 偶数 时 ， fal ae 而 1 为 奇数 时 ， tt 

由 于 我 们 考察 的 问题 是 分 伟 点 附近 极 大 值 线 的 拓扑 情况 ， 所 以 对 s, t 做 一 些 限 制 
是 必要 的 。 

令 ssOg)h t-O(e) db — 0 0000 eee (6.4.2) 

TE LABES RISE, TSAA FEES HÉT AT 两 个 次 数 不 超过 1 一 1 阶 的 多 项 陈 
EA "IBU 函数 。 


locu dius c des] "m 
Ty - 4", アニ ーー, [Ham eee 6.4.3 
a 5» 2i } 2 人 ) 
| LIT titl (o l- 
T. = -4p Ne get l = 700703 LAR, iau. 6.4.4 
6 5x. 2i 1 | 2 (6.4.4) 


讨论 ， 当 ;7 为 偶数 时 ， (644) 式 中 t= 2^ BURT, = 0 的 一 个 解 。 在 这 种 情况 下 ， 
Ei Aa > 0 或 Aa <0， 一 条 极 大 值 线 都 会 沿 着 两 个 方向 发 展 ， 而 当 ! 为 奇数 时 ， 由 于 
1 始终 为 偶数 次 守 ， 在 Aa <0 时 ，T 就 会 减 小 。 在 这 种 情形 下 ， 随 着 a 的 增 大 ， 分 贫 
点 处 的 极 大 值 线 也 就 不 可 能 产生 分 裂 。 

在 Aa <0 时 ， 为 了 得 到 分 岔 点 附近 的 近似 解 ， 根 据 (6.4.2) 式 的 要 求 ， 不 妨 令 
= Cr ， 代 入 (6.4.3) 式 和 (6.4.4) 式 得 到 


| Ld ied yu 3 m 
T. ecc bpm > p = 0, ] ET 
l a ML 1 2i | 2 


1 NE du 
T, = php n 7 es ci Ns 7 こさ 
^ à m per | 


r=0 i] 


对 于 每 一 个 1，C 不 依赖 于 2 和 a， 给 C 赋 值 ， 就 可 以 确定 一 条 抛物 线 。 在 平面 
(s,?) 内 ， 在 一 个 分 合 点 (a6,5,) 处 交汇 的 所 有 极 大 值 线 都 是 抛物 线 ， 由 曲率 和 点 值 就 
可 以 确定 这 些 极 大 值 线 在 (a,,b) 附 近 的 分 布 情 况 。 


6.5 基于 小 波 变换 模 极 大 值 的 多 重 分 形 谱 的 计算 


从 以 上 分 析 可 以 看 出 ， 小 波 变换 模 极 大 值 线 (已 结构 已 经 包含 了 奇异 性 的 必要 信 
息 ， 但 由 于 分 形 的 奇异 点 是 非 孤 立 性 的 ， 要 计算 点 态 奇 异 指数 cc(xo) 是 很 困难 的 。 为 
此 ， 我 们 借助 由 Ameodo、Bacry 等 人 89 引入 的 一 个 和 函数 S(g,a) (该 函数 可 被 用 来 
测量 所 有 小 波 变换 模 值 q 次 窜 的 和 ) 通过 小 波 变 换 局 部 极 大 值 来 测量 多 分 形 谱 。 


ts 分 形 理论 中 若干 问题 的 小 波 解 法 


S(a, a) = Y, ua, b, y 
íeL, 


RPW, ula, b j~ a^ ， 玉 代表 在 尺度 a 下 第 ;个 极 大 值 线 的 位 置 。 
计算 多 分 形 谱 万 的 经 典 方法 是 由 下 式 估算 r(g): 
S(g,a) ~ a" 
然后 通过 Legender 变换 得 到 -OY _1) 。 然 而 在 实际 应 用 中 由 于 数据 集 噪声 和 
有 限 长 度 的 缘故 ， 就 会 有 许多 问题 产生 。 为 此 我 们 将 做 下 面 的 两 项 工作 。 

在 每 一 尺度 a 下， 对 于 某 个 a > 0， 人 和 假设 任何 相 邻 的 极 大 值 点 b 、b,,, 有 一 段 距 
Bib, -Bb|>s， 否 则 在 区 间 e 上 , S(a,qg) 仅 包含 最 大 幅度 的 极 大 值 ， 这 样 S(a,g) 就 
不 至 于 由 于 快速 震荡 产生 的 较 近 极 大 值 引发 的 数值 又 增 。 

对 于 S(q,q)， 当 g < 0 时 ， 由 于 极 大 值 线 上 点 附近 的 小 波 变换 模 值 很 小 ， 而 
Sa,9) 值 主要 依赖 于 较 小 幅度 的 小 波 变换 模 值 ， 这 样 S(a,q) 计 算 就 缺乏 稳定 性 。 为 
了 避免 这 种 不 稳定 性 ， 对 于 小 于 a 的 所 有 尺度 上 每 一 个 极 大 值 点 处 的 小 波 变换 模 极 
大 值 用 沿 对 应 极 大 值 线 小 波 变换 模 极 大 值 的 上 确 界 代替 。 

在 每 一 尺度 a 下 ， 我 们 定义 一 个 类 似 于 “概率 ”的 量 P(g,a) 


Pl(q,a)= Pe QA a) 


于 是 广义 次 信息 量 维 数 Dla) A: 
1 YP (aa)inP(a.a) 


D(q) = lim 一 一 一 一 一 一 ， (g #1) 
a30 gq —] Ina 


再 由 第 二 节 有 关 推 导 结 果 ， 求 得 奇异 性 指数 a(q)= 一 [(q -1}Dp(q， 同 时 多 分 形 
AREN (分 形 集 的 Hausdorff HER) A: 


f(a)= qa(q)- (a - の (2) 


d 
dq 


6.6 仿真 算 例 分 析 


例 1: 在 小 波 变换 极 大 值 处 理 方面 ， 分 形 信号 的 奇异 性 和 一 般 信号 奇异 性 有 所 不 
同 。 图 6.1(a) 为 原始 波形 信和 号， 在 [50.100] 一 段 包 含 了 一 个 分 形 结构 ， 采 用 高 斯 函数 
作为 母 小 波 ， 对 其 进行 小 波 变 换 极 大 值 分 析 。 图 6.1(b) 显 示 了 两 类 信号 不 同 的 表现 : 
在 1=160 处 的 小 波 变 换 极 值 在 各 尺度 下 都 可 以 很 好 地 跟踪 ， 而 在 := 70 处 的 小 波 变换 
PHATE j < 4 的 情况 下 分 裂 成 多 个 极 值 点 ， 跟 踪 问 题 由 此 产生 了 ，。 
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例 2， 以 概率 p =0.4 和 っ ヵ , = 0.6 的 魔鬼 阶梯 《 见 图 6.2(a)) 为 例 ， 用 文中 方法 
计算 其 谱 维 数 。 小 波 函 数 y = -9 ，9 为 高 斯 函数 。 图 6.2(b) 显 示 基 于 极 大 值 方法 的 
多 分 形 谱 (+ 号 曲线 ) f(g) 和 用 Legendre 变换 方法 所 计算 的 谱 比较 接近 。 


6.7 本 章 小 结 


在 本 章 ， 首 先 介绍 了 多 分 形 谱 的 两 套 描述 语言 : w- f(g) 语言 和 g -DD, 语 言 ， 分 
析 了 现今 有 关 多 分 形 谱 尊 异性 刻画 方面 存在 的 问题 ， 提 出 了 基于 高 斯 小 波 函 数 跟踪 
极 大 值 线 的 一 种 尺度 参数 化 (a( 站 2 吉方 法 。 

在 第 四 节 里 ， 针 对 极 大 值 分 岔 点 附近 的 性 态 ， 引 入 了 两 个 多 项 式 函 数 荆 、 卫 作 
AF Ra, 探究 了 出 現 分 盆 点 的 可能 | 

文中 第 五 节 ， 多 分 形 谱 奇异 性 的 几乎 所 有 信息 都 包含 在 极 大 值 线 结构 中 ， 为 此 
我 们 仿照 信息 量 维 数 的 定义 ， 基 于 小 波 变换 模 极 大 值 ， 提 出 了 多 分 形 谱 的 一 种 估算 
方 法 。 

最 后 ， 列 举 了 两 个 范例 。 其 中 第 一 个 例子 采用 了 包含 一 段 分 形 结构 的 信号 ， 经 
小 波 极 大 值 处 理 后 ， 从 图 上 可 以 看 出 ， 分 形 的 极 大 值 点 出 现 分 裂 ， 有 必要 进行 跟 
踪 。 第 二 范例 则 以 一 个 魔鬼 的 阶梯 为 例 ， 显 示 了 文中 提出 的 谱 估 算 方 法 是 有 效 的 。 


在 多 分 形 谱 的 研究 过 程 中 ， 本 人 认为 ， 尚 有 必要 在 以 下 方面 进一步 深入 下 去 : 
要 正确 地 揭示 多 分 形体 时 空 演变 的 动力 学 过 程 ， 和 本 论文 第 四 章 提 出 的 单 分 形 
时 变 维 数 一 样 ， 有 必要 对 多 分 形 引 入 时 变 维 谱 的 概念 ， 可 考虑 从 下 面 的 两 个 途径 开 
展 研究 : 
C1) 仿照 现代 信号 处 理 中 有 关 时 变 谱 的 分 析 理 论 ， 用 小 波 分 析 手 段 处 理 时 变 维 
i. 


(2) 采用 Marcov 过 程 理论 对 多 分 形 时 变 维 谱 的 漂移 加 以 刻画 。 


分 形 理论 中 者 干 问题 的 小 波 解 法 


r(t) 
i 
0 50 100 150 200 250 


图 6.1(a) 包含 [50, 100] 一 段 分 形 信号 的 原始 波形 


WT, x (1 )-À — — Ate $< gt par 
WT;? x (2) タス ーー ペー ター ター 
WT,3 ri 会 一 -一 一 一 一 -一作 < 一 一 一 一 一 
WT, ri 一 -一 一 一 一 一 -一 一 
WT。s rz) 

gat eee 


) 
0 50 


图 6.1(b) [50，100] 一 段 难以 进行 小 波 变换 极 大 值 跟踪 


z 
. "EN 
1000 2000 3000 4000 


6.2(a) p, = 0.43 p, = 0.6 的 魔鬼 阶梯 


0.6 oe 
E 
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图 6.2(b) 基于 极 大 值 方法 的 多 分 形 谱 /(〆 ) 和 用 Legendre 方法 
所 计算 的 谱 ， 前 者 为 + 号 曲线 ， 后 者 为 实 曲线 。 
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小 波 理论 中 的 多 分 辨 分 析 思 想 ， 体 现 了 人 们 认识 和 识别 形体 的 过 程 遵循 了 一 
种 从 低 分 辩 率 到 高 分 辩 率 的 原理 。 对 于 分 形 的 观察 也 正 是 这 样 的 方式 。 从 大 到 小 
的 尺度 变换 ， 在 越 来 越 小 的 尺码 上 观察 越 来 越 丰富 的 分 形 细节 问题 。 小 小 分 析 和 
分 形 具 有 深刻 的 内 在 联系 ， 两 者 在 尺度 变换 上 具有 一 致 性 ， 分 形 是 一 种 几何 语言 ， 
而 小 波 是 一 种 分 析 局 域 奇异 性 和 非 平稳 性 的 工具 ， 两 者 相得益彰 。 

小 波 从 诞生 到 现在 不 过 十 几 年 的 时 间 ， 已 经 有 了 坚实 的 理论 基础 和 广阔 的 应 
用 领域 ， 而 分 形 理论 和 小 波 相 比 虽然 起 步 比较 早 ， 但 因 其 “ 离 经 叛 道 ”的 反日 性 
和 独立 性 不 易 为 大 多 数 的 人 们 所 接受 ， 况 且 本 身 的 理论 基础 尚 建立 得 不 够 牢固 ， 
发 展 研究 因 些 就 相对 显得 较为 缓慢 。 小 波 分 析 工 具 从 开始 进入 分 形 领域 至 今 则 短 
短 不 到 十 年 的 时 间 ， 取 得 了 一 些 可 喜 的 成 果 ， 但 分 形 理论 中 有 许多 存在 的 问题 ， 
还 有 待 借助 包括 小 波 分 析 在 内 的 方法 去 解决 ， 已 有 基于 小 波 方法 的 结果 也 有 个 足 
之 处 。 但 我 们 可 以 确信 的 是 ， 小 波 分 析 的 引入 犹如 催化 剂 增添 了 分 形 研 究 的 活力 ， 
能 够 为 分 形 理 论 的 进一步 发 展 和 应 用 起 到 一 定 的 促进 作用 。 

基于 以 上 理念 和 思路 ， 在 前 人 工作 的 基础 上 ， 本 文 利 用 小 波 工 具 就 分 形 理论 
有 关 问 题 进行 了 较为 深入 的 探讨 。 

本 文 的 主要 工作 和 创新 之 处 可 归纳 为 以 下 儿 个 方面: 

1. 分 析 了 分 形 曲线 维 数 的 传统 计算 方法 及 其 缺陷 ， 提 出 多 分 维 分 形 曲 线 的 
模式 。 

2. 基于 小 波 阔 值 去 噪 方法 ， 给 出 曲线 上 点 的 局 部 分 形 维 数 的 计算 公式 及 算 
法 。 通 过 收敛 性 的 证 明 ， 该 方法 可 以 有 效 地 保留 其 分 形 维 数 。 仿 真 计算 
结果 表明 该 方法 是 有 效 的 。 

3. 在 处 理 分 形 反 问题 方面 ， 和 前 人 所 选 特征 点 不 同 ， 我 们 的 做 法 是 选择 小 
波 变换 零点 作为 特征 点 ， 处 理 对 象 定位 于 分 形 插值 函数 。 

4. 基于 小 波 变换 在 零点 具有 的 的 保 相 似 形 特征 ， 在 已 知 分 形 插值 函数 图 象 
情况 下 ， 从 有 限 数目 个 零点 出 发 ， 给 出 了 寻找 一 组 仿 射 压缩 变换 使 该 冰 
数 图 象 成 为 未 知 和 迭代 函数 系统 吸引 子 的 一 种 构造 方法 。 算 例 结 采 也 验证 
该 算法 的 有 效 性 。 

s. 提出 了 分 形 时 变 维 数 的 概念 ， 并 将 其 引入 两 类 具有 局 部 自 相 似 性 的 有 侦 
随机 过 程 中 。 

6、 借 助 小 波 基 固 适合 于 分 析 非 平稳 性 和 局 部 自 相似 性 这 一 特点 ， 给 出 Hurst 
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相 容 的 ， 并 用 仿真 结果 加 以 验证 。 

7， 利 用 时 变 Hurst 指数 理论 对 上 海 股票 市 场 周 指数 变化 进行 了 实证 分 析 ， 
结果 表明 该 理论 对 投资 决策 有 指导 意义 。 

8， 分 析 了 多 过 程 的 小 波 变 换 系数 之 间 在 相关 结构 上 的 特点 以 及 小 流 消 失 
矩 的 大 小 对 系数 相关 性 衰减 性 的 影 啊 。 

9， 借 助 小 波 变 换 对 自 相似 过 程 所 起 的 白化 作用 和 最 小 二 乘法 ， 在 无 须知 
道 yy 过 程 中 参数 分 布 的 情况 下 ， 给 出 了 基于 小 波 系数 方差 结构 的 谱 参 
数 的 估算 式 。 用 仿真 对 方法 的 有 效 性 进行 了 验证 ， 结 果 较 为 满意 。 


10， 提 出 在 最 优 尺 度 上 跟踪 小 波 模 极 大 值 线 的 参数 方法 ， 为 研究 多 分 形 庶 
奇异 性 提供 了 一 种 途径 ， 可 节省 大 量 的 运算 量 。 

11. 选取 了 小 波 变换 在 分 岔 点 邻 域 Taylor 级 数 展开 中 的 两 个 多 项 式 函 数 作 
芽 函 数 ， 探 讨 了 分 形 测 度 小 波 变换 极 大 值 线 的 拓扑 分 贫 情 况 。 

12， 提 出 了 一 种 基于 小 波 变 换 模 极 大 值 的 多 重 分 形 谱 估 算 方法 ， 以 “魔鬼 
的 阶梯 ”为 例 对 该 方法 的 有 效 性 加 以 验证 ， 取 得 了 和 Legendre 变换 万 


法 几乎 相同 的 结果 。 
在 分 形 理论 的 研究 过 程 中 ， 本 人 认为 ， 有 必要 在 以 下 儿 个 方面 沿街 进 一 步 深 


入 探讨 : 

|， 为 了 能 够 完整 地 刻画 多 分 形体 在 时 空中 演变 的 动力 学 过 程 ， 有 必要 引入 
时 变 维 谱 的 概念 。 仿 照 现 代 信 和 吴 处 理 中 有 关 时 变 谱 的 分 析 理 论 ， 用 小 六 
分 析 手 段 处 理 时 变 维 谱 。 

2. 采用 Marcov 过 程 理论 对 多 分 形 时 变 维 谱 的 漂移 加 以 刻画 。 

3， 由 小 波 变换 值 的 关系 可 以 反 推 分 形 的 变换 关系 ， 然 而 ， 这 种 反 推 是 否 确 
定 唯一 的 分 形 ? 这 是 一 个 尚 待 解决 的 问题 。 

4， 对 于 非 相 似 分 形 《〈 如 圆 和 椭圆 、 球 和 椭 球 ) 的 小 波 分 析 反 演 ， 这 必然 水 
及 尺度 参量 应 该 为 n 维 对 角 和 矩阵 的 小 波 函 数 。 

g(x)= Y'ap(Ax- b), AAT BER. 
bez? 


有 关 这 种 小 波 函 数 的 研究 已 有 不 少 。 尽 快 地 将 其 用 于 这 类 分 形 ， 这 也 
是 以 后 反 演 工作 的 一 个 研究 方 网 。 
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